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FOREWORD 

I 
! 

j The program documentation presented herein is a follow-on effort to 

that described in "136 MHz/400 MHz Earth Station Antenna-Noise Temperature! 
Prediction Program for RAE-B", by Ralph E. Taylor, Joseph J. Fee and M. 

t 

Chin, NASA/GSFC Report No. X-752-72-324, September 1972. Sections 1.0 and 

| — - 

2.1 are taken directly from this reference. 
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SECTION 1.0 
INTRODUCTION 


In 1973, the Radio Astronomy Explorer-B (RAE-B) 
satellite will be placed in a 1100 km-altitude circular 
orbit around the Moon to make radio astronomy measurements. 

The purpose of the simulation study described in this 
report is to determine the 136 MHz and 400 MHz noise temper- 
ature of the ground network antennas which will track the 
RAE-B satellite during data transmission periods. Since the 
noise temperature of the antenna effectively sets the Signal- 
to-Noise Ratio (SNR) of the received signal, a knowledge of 
SNR will be helpful in locating the optimum time windows for 
data transmission during low-noise periods. 

Antenna-noise temperatures at 136 MHz and 400 MHz will 
be predicted for selected earth-based ground stations which 
will support RAE-B. Telemetry data acquisition will be at 
400 MHz; tracking support at 136 MHz will be provided by the 
Goddard Range and Range Rate (RARR) stations. 

The antenna-noise temperature predictions will include 
the effects of galactic-brightness temperature, the sun, and 
the brightest radio stars. Predictions will cover the ten- 
month period from March 1, 1973 to December 31, 1973. The 
RAE-B mission will be especially susceptible to SNR degrad- 
ation during the two eclipses of the Sun occurring in this 
period. . • , 

The RAE-B Tracking Antenna Noise Temperature Program 
Documentation is to define the operation of the prediction 
program in terms of the mathematical operations, input data 

requirements and capabilities of. the program. Many, of the 

| 


Program Modules utilized' were previously developed for the 
Data Quality Prediction Program (DQP) developed under NAS 
5-11736 MOD 106 and if unchanged, their descriptions are 
referenced in that program rather than repeated here. Each 
new module or subroutine and any changed routines are fully 
documented in Section 3.0 of this report and a flow chart 
of the overall program is presented. 



SECTION 2.0 
PROGRAM DESCRIPTION 

• • • , . .... ' .. r 

2.1 DEVELOPMENT OF EQUATIONS 

Four sources of antenna -noise temperature are con- 
sidered in this study: 


o • Sky-Brightness Temperature 

o Sun 

e Radio Stars 

• o Antenna Back Lobe Noise Temperature 

e Total Antenna-Noise Temperature 

The formulation utilized within this program has 
been previously utilized in the Data Quality Prediction 
Program [1] developed by Wolf Research § Development Corp. 
for NASA/GSFC under contract NAS 5-11736 DCN 523-W- 70446. 
The equations presented in the following are taken from 
References 1 and 2. 

a) Sky-Brightness Temperature 


Kraus [3] develops the following formulation for 
sky brightness temperature 


G“90° - 0 o 4*“ 2 Tr 

f j T (0,4>) G (6,4)) sin 9303<I> 

Jq 


T SKY- 


0 = 90° - 8' 0 4^“ 2ir 

f / G (6,4)) sin 63034) 

J o Jo 


3 


where 


®0 = elevation angle between antenna's boresight 

axis and the horizon, degrees 

T (0 ,<{>) is the noise temperature distribution o£ the 
galaxy (excluding the sun and predominant 
radio stars) obtained from References 1. and 2. 

G(6,4>) is the lossless antenna gain distribution. 

0 and <}> are the orientation angles defining the 

position of a radial surface element within 
the celestial hemisphere. • 

Figure 1 shows the physical relationshiD of the 
variables given in the above equations. The double integral 
is computed by rectangular integration assuming- that the 
antenna boresight is directed at the center of the Moon's 
optical disk. The antenna patterns and brightness tempera- 
tures are accessed from magnetic tape storage in the computer 

program. 


4 




Pqlt-ern to Celestial Sphere 


figure l 


Relation of Anlcnna 



This program required the use of an accurate radio- sky 
map which covers the celestial sphere completely, for both 
136 MHz and 400 MHz. Since detailed sky-brightness temperature 
contours were not available at these frequencies, it became 
necessary to scale either existing radio maps in temperature, 
or to generate a composite map from various smaller maps. The 
136 MHz and 400 MHz sky-brightness temperature maps , were pre- 
pared in this manner. 


The 136 MHz radio map was scaled from data at 150 MHz, 
published in 1971 by Landecker and Wieleb inski (Reference 4). 
The following relationship was used for scaling. 



Reference 5 was utilized for the conversion of galactic 
coordinates, employed by Landecker and Wieleb inski, into the 
necessary equatorial coordinates required for this program. 

The 400 MHz radio map is a composite map formulated 
frem the sectional maps published in 1962- by Pauliny-Toth 
and Shakeshaft (Reference 6) , and in 1956 by Droge and Priester 

(Reference 7) . 


b) Sun 

The contribution of the sun to the antenna-noise 
temperature is given by Berkowitz [8] as: , 


T SUN 



T,, assuming 6 >> 0 

b ’ a s 


6 



where 


6 g = Angular radio diameter of sun^ apparent temper 
ature model, degrees; assume & s = 0.66° at 1-36 
MHz and 400 MHz. 

9 = Half-power beam width (HPBW) of symmetrical 

antenna main lobe. 

= 8 x 10^ K for quiet sun ideal model at 136 MHz, 

and 6 x 10 5 K at 400 MHz. 


c) Radio Stars . . 

The following equation is used by Taylor '[1,9] to 
compute antenna noise power rise due to a point-source radio 
star: 


N, 


M 2 

1 G G (0) a 

_ _P D Af 

2 4ir 

n=l 



for a single polarization. 


where 


M - is the number of radio stars 



is the observed radio star noise flux density, 
7 - 1 

wm Hz , constant over bandwidth Af. 
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X - wave length of transmission 

Gp - peak antenna power gain, above isotropic 

G(8)- antenna gain attenuation. at angle 6 of f-boresight 
i.e. , G (6) = 1 for 0 =0 . 

Note that the antenna-noise temperature is 

iT*- - 

N* 

T* = degs.. K 

kAf 

where 

N* = total noise power due to all radio stars- 
within the antenna’s radiation pattern 

' -23 

k = Boltzmann’s constant, 1.38x10 J/K 
Af = noise bandwidth of receiver, Hz 

d) Antenna Back Lobe Temperature 

The black-body radiation of the Earth contributes to 

the overall antenna-noise temperature by means of the back 

* * 

lobe of the ground antenna. 

Based on Blake’s data (Reference 10), the effect of 
antenna back lobe temperature, T BACK , is approximated by 
adding a constant to the equation for antenna-noise temper- 
ature as follows: 


• 8 




136 

MHz 

t back ~ 

75°K 

400 

MHz 

T BACK = 

35°K 


e) Total Antenna-Noise Temperature 

The total antenna-noise temperature, T TQT , is computed 
by summing each of the four contributions. 

T TOT = t sky * t sun + T * + t back 

Antenna-noise temperature will be maximum at New Moon, 
once each month. A higher peak is reached once each year 
(December) when the Galactic Nucleus is eclipsed by the Sun, 
during New Moon (see Figure 2). 
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(New Moon To 
New- Moon) 
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Earth 


\ 


\ 


New 
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(Position: December 24, 1973 
e 15 h 07” ZULU) 
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Figure 2. Celestial Source Spatial 
Arrangement in Month of December 



2.2 PROGRAM LOGIC FLOW 

The following logic steps determine the antenna-nois 
temperature for each ground antenna which will track the 
RAE-B: 

e Eafch station is tested for a visible moon. 

® Antenna-noise temperature is computed for 
each. station having a visible moon. 

© Antenna-noise temperature is recomputed at 
periodic time intervals from Moonrise to 
Moonset . 

The logic of the program will be described by an 
overall program flowchart presented in the following: 
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SET DATE FOR OR3IT CALCULATIONS 
ANDSET JSECN TO ZERO TO KEEP 
TRACK OF PREDICTION TIME INTERVAL 



CHECK FOR STOP TIME 


CHECK FOR ONE WEEK PASSBY 
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YES 



IF IRSINE IS NEGATIVE, MOON IS BELOW 
HORIZON 


COMPUTE STATION-MOON UNIT VECTOR 
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ISTP = ISECN.GT.ITIME 




RESET NEWSTART DATE FOR NEXT 
PREDICTION PERIOD 














•SECTION 3.0 

PROGRAM DESCRIPTION AND LISTINGS 


3.1 PROGRAM REFERENCE CHART ■ 

Name in (In DQP 

Name the DQP Page Documentation) 


MAIN PROGRAM (ATEMP) 

DQP 

> 

1 

* 

BLOCK DATA 

BLOCK DATA #1 

A- 31 

ADDYMD 

ADDYMD 

A- 4 7 

ANGLE 

ANGLE 

A- 51. 

ARCTAN 

ARCTAN 

A- 54 

BASEYR 

BASEYR 

A- 57 

CROSS V 

CROSSV 

A- 90 

DIFF 

DIFF 

A-92 

DIFFTM 



DINRAD 

DINRAD 

A-97 

DOTTED 

• DOTTED 

A-100 

EPHEM 



HEADIN 



MODULO 

MODULO 

A-lll 

ORBFIX 

ORBFIX 

A-130 

ORBMS 



OUTPUT 



PRTOUT 



REVERS 

REVERS 

A-187 

ROTATE 

ROTATE 

A-190 

RYMDI 

RYMDI 

A-193 

SOURCE 

NO I SPIV 

A-114 

STATIO 

STATIO 

A-207 


TP LOT 
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ROUTINES WITH CHANGES FROM DQP 


Name 

Page 

Name in 
the DQP 

Page 

ATEMP 

27 

DQP 

A - 1 

BLOCK DATA 

46 

BLOCK DATA n 1 

A- 31 

SOURCE 

48 

NOISPW 

A- 114 


NEW ROUTINES 


Name 

Page 

EPHEM 

55 

HEADIN 

59 

OUTPUT 

62 

PRTOUT 

64 

TPLOT 

67 
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3.2 COM’ 10 M BLOCKS CROSS RHFLRiLMCl: CHARI 


COMMON BLOCKS 



2 

E 

■ _j 
< 

c 

2 

1- 

2 

< 

_j 

ea 

> 

< 

□ 

INERTL j 

MISCEL j 

NEW J 

j 1N3IU0 

RECEV I 

ROTATS 

STAID j 

ATEMP 

■ 


B 

B 

B 

B 

B 

B 

B 

B 

ADDYMD 



B 





IB 


B 

BLOCK DATA 

9 

B 

B 





B 



DIFF 



B 








DIFFTM 



B 

1 


| 


Hi 


|| 

HEADIN 

m 

B 

B 



B 


B 

| 

B 

ORBFIX 

u 



B 



1 




ORBMS 



B 

B 

n 


B 

Hi 


H 

j PRTOUT 

m 




B 



Hi 

Brail 

HI 

SOURCE 





B 

B 

X 

B 

B 

B 

ST ATI 0 





B 

B 

B 

B 

B 

B 

TPLOT 

IB 

IB 

B 



B 




B 
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DESCRIPTIONS OF ROUTINES WITH CHANGES FROM DQP 

IDENTIFICATION 


A. Name: ATEMP 


B. Language: FORTRAN IV 


C. Machine: IBM 360 


D. Purpose: To predict the RAE-B tracking station 

antenna-noise temperature at both 136 
megacycle frequency and 400 megacycle 
frequency for one week period as well 
as ten months period. 


E. Common Blocks : 


Name Variable Dimension Type Description 

/INERTL/ ORBELM (6) R*8 Input-Inertial rec- 

tangular coordinates 
(1-3) and velocities 
(4-6) of satellite 
in meters and meters 
per second respectively 

THDOT (1) R*8 Output-Rotational 

rate of earth in 
radians per second 

/RECEV/ INTGTH (1) R*4 Input-The orientation 

angles increment of 
6 (in deg. ) 

INTGPH (1) R*4 Input-The orientation 

angles increment of 
(in deg. ) 

MESHRC (18732) 1*2 Output- Containing 

receiver antenna 
gains and sky map 

/ROTATS/ I TYPE (1) 1*4 Output-selected track- 

ing system type code 
(1=SATAN system, 2= 
85-foot dish, 5=40- foot 
- dish, 4=16- element yagi 
array, 5=19- element 
yagi array) 
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Name Variable 

Dimension 

T YP.£ 

Description 

MATRX 

(21) 

1*4 

Temporary locations 

/STAID/ NAME 

(7) 

OO 

•K 

Output- Tracking station 
identification 

TYPE 

(7) 

1*4 

Output-Tracking system type 
code 

LATD 

(7) 

1*4 

Output-Degrees of station 
geodetic latitude. 

LATM 

(7) 

1*4 

Output -Minute's of station 
geodetic latitude 

LSLAT 

(7) 

R*4 

Output- Seconds of station 
geodetic latitude 

LOND 

(7) 

1*4 

Output - Degrees of station 
geodetic longitude 

LONM 

(7) 

1*4 

Output-Minutes of station 
geodetic longitude 

LSLON 

(7) 

R*4 

Output- Seconds of station 
geodetic longitude 

LGAINA 

(7) 

R*4 

Output-Receiver's peak antenna 
gain (dB) 

LELMIN 

(D 

R*4 

Output-Minimum observational 
altitude angle (deg.) 

NSTA 

(1) 

1*4 

Output-Number of tracking 
stations 

LEIGHT 

(7) 

R*4 

Output -Topographic height 
of stations (meters) 

/ORIENT/ LHATS 

(42) 

R*4 

Input-Local north, east vectors 
at. tracking sites 

LZHAT 

(3,7) 

R*4 

Input-Station's local zenith 
unit vector 

I SUN 

\ 

(2) 

1*4 

Temporary location 

1 

I WAVE 

(D 

R* 4 

Output- Selected transmission 
wave length (meter) 
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Name 

Variable 

Dimension 

Type 

Description 


I BAND 

CD 

R*4 

Input-Selected transmission 
band width (H^) 

/MISCEL/ 

EPOCHY 

CD 

1*4 

Output-Date (year , month , day) 
o£ orbit's epoch (YYMMDD) 


EPOCHH 

CD 

1*4 

Output-Time (hour , minute , 
-second) of orbit's epoch 
(HHMMSS) 


INTRVL 

CD 

1*4 

Output- Time interval of orbit 
points (seconds) 


ISTRTY 

CD 

1*4 

Output-Date of start of pre- 
dictions (YYMMDD) 


ISTRTH 

CD 

1*4 

Output-Time of start of pre- 
dictions (HHMMSS) 


NUMSAT 

CD 

1*4 

Output-Satellite sequence 
number 


LTSATX 

CD 

R*4 1 



LTSATY 

CD 

R*4 

> Input-Station- to-satellite 


LTSATZ 

CD 

R*4 ] 

| vector components (meters) 


LSTAX 

CD 

R*4 ^ 

1 


LSTAY 

LSTAZ 

CD 

CD 

R*4 
R*4 j 

j> Input-Station's geocentric 
coordinates (meters) 


TA 

CD 

R*4 

Input-Predicting temperature 


IRANGE 

CD 

R*4 

Output -St at ion- to-satellite 
distance (meters) 


IRSINE 

CD 

R*4 

Output -Sine of satellite 
elevation angle 

/NEW/ 

FREQ 

F. Name 

CD 

List : 

1*4 

Output- Frequency (MH^) 

Name 

Variable 

Dimension 

Type 

: Description 

/CONTL/ 

TAPE21 

CD 

L 

Output-Code indicating the 
presence of satellite transmi 
antenna gains input tape on 


unit #21 
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Name Variable 

Dimension 

Type 

Description 

/INTGER/ INTGTH 

(1) 

R*4 

Output - angular 

increment 




of 8. (deg.) 


INTGPH 

(1) 

R*4 

Output- angular 

increment 

** 

*> 



of <t>. (deg.) 



G. Non-System' Routines Required: 


PLOTST, CALCOM, DIFF, BASEYR, STATIO, HEADIN, ROTATE, 
ORBMS , ORBFIX, ADDYMD , REVERS, DOTTED, PRTOUT 


H. Input/Output Units Required: 


Logical Unit 

Type 


Data 

FT04F001 

Disk Input 

Ephemeris routine 

FT05F001 

System 

Input 

Job control instructions, 
input data sets 

FT06F001 

System 

Output 

Printed data sets 

FT23F001 

2400-9 

tape 

Cumulative data 

PLOT APE 

2400-9 

tape 

Output for the Cal comp 
Plotter 

FT21F001 

2400-9 

tape 

Spacecraft Antenna gains 


and- brightness sky temperature 
data. 

II. METHOD 

This program is divided into four major operational 
sections : 

1. Access and storing of control information, 
and tracking systems parameters. 

2. Access and storing of moon parameters. 

• 3. Moon orbit, calculations and storing. 

4. Telemetry data predictions and display. 
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d. Temperature for quiet sun ideal model 

e. Angular diameter of sun's apparent temperature model 

f. Radio stars identification 

2. Predicted temperature data 

* 

a. Time of prediction 

b. Predicted temperature 

c. Tracking station 

d. Temperature source 

The entire operational sequence described in this section is 
repeated for each tracking station, when this is completed, 
the process return to section 3 in search for the next tracking 
period. 


XII. 

CONSTANTS AND i 
A. Constants 

MESSAGES 


Name 

Value 

Dimension 

Description 

SPACE 

10080 

CD 

Storage space reserved for 
orbit computations 

AUNIT 

6378166. 

6378388. 

6378165. 

(3) 

Equatorial radii of earth 
ellipsoids in meters 

EUNIT 

298.25 

297.00 

298.30 

(3) 

Reciprocals of flattening 
of earth ellipsoids ■' 

GMUNIT 

.0012394454711 
. 0012394153954 
.001239443589 

(3) 

Square roots of the products 
of the earth's masses and 
gravitational constants in 
units of (earth radii) exp 
3/2 per second 
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Name 

Value Dimension 

Description 

LIGHT 

2.997925X10 8 

CD 

Speed of light in meters/ 
second 

RD 

.0174532925 

CD 

One degree expressed in 
radians 

TWOPI 

6.28318530-72 

CD 

2 T 

THDOT 

. 729211585648X10' 

4 (D 

Rotation rate of earth in 
radian/second 


B. Messages: 




Message 


Description 


INVALID INPUT CODE • This message is printed when the 

tracking system type code, satel- 
lite type code, and earth constants 
code could not be recognized. The 
data set is bypassed. 

INPUT TAPE TRANSMISSION When the value of control para- 

ERROR meter TAPE 21 is .TRUE, and read 

attempts on logical unit 21 are 
■unsuccessful, this meeting is 
printed and the process is 
terminated . 

RECEIVER ANTENNA GAINS The message is self-explanatory. 

TABLE IS NOT AVAILABLE Noise power will not be predicted. 

FOR STATION XXXXXX - 
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Section 1 and Section 2 are performed only once, while 
program execution returns to the other sections repeatedly for 
temperature predictions from various prediction period. In 
the following the logical structure of each program section 
will be discussed. 


Section 1 


The major functions of Section 1 are to: 

a. Access and process internally stored control 
information (select and initialize display 
hardware, set up constants). 

b. Access and process inputted control data 
(initialize input unit of transmitter antenna 
gains, store prediction time). 

■. 'i 

c. Access, store, and display tracking system 
parameters which are: 

1. Station identification (6 alphanumeric 
characters) 

2. Tracking system type codes: 

* 

1 . = SATAN sys tern 

2 = 85 -foot dish 

3 = 40-foot dish 

4 = 16 element Yagi array 

5 = 9 element Yagi array 

3. Station's position [geodetic latitude, 
longitude (degree, minute, second), 
topographic height (meters)] 
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4. Minimum look angle above horizon (deg.) 

5. Receiving antenna gain above isotropic 
source (dB) 

► 

Section 2 

Moon parameters consist of : 

1. Identification (6 alphanumeric characters) 

2. Keplerian orbital elements: 

a. Epoch of elements (Date: year, month, 

day; Time: hours, minutes, seconds) 

3. Earth constants code: 

0: a = 6378166m, f=l/298.25 

1: a = 6378388m, f=l/297.0 

2: a e = 6378165m, f=l/298.3 • 

4. Transmitter frequency (megacycle/second) 

To assure consistency between orbital elements and 
station coordinates, tracking station positions are computed 
on an ellipsoid model specified by the earth constant code. 

Section 3 

Substantial time can be saved in computing efforts, by 
calculating and storing moon ephemerides for up to one week 
predictions. Using the JPL ephemeris routine, the inertial 
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orbital positions and velocities are obtained analytically 
(subprogram ORB'IS) , which are then converted to a geocentric, 
Greenwich oriented reference system (subprogram ORBFIX) . This 
sequence is repeated at each incremented time point returning 
to the orbital calculations section of the process. Computer 
memory is reserved for the storing of orbits for the period 
of one week at 30 minute intervals. 


Section 4 


Telemetry data is predicted for each tracking site from 
a given antenna. Each stored orbital point is tested for 
visibility before utilizing it in the prediction process by 
forming the scalar product of station's zenith vector and station- 
to-moon vector. When the orbital point is above the minimum 
look-angle, slant range and transmitter antenna gain are in- 
cluded in the predicted temperature. The temperature data is 
obtained by calling subprogram SOURCE, and its value reflects 
the effect of galactic, solar and major radio stars. The 
position of the sun is considered constant within the time of 
moon passage consequently, the right ascension and declination 
of sun are computed at one hour intervals. 

Predicted temperature data in display on the computer 
system output unit in printed form, when the predicted time 
is equal or more than one week it is also plotted by using 
WOLF calcomp plot package. The following items comprise the 
printed output: - 

1. Identif ication heading 

a. Tracking period and interval 

b. Transmitter frequency (MH^) 

c. Station identification 
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o 

o 

2 > 

t 

.3 

0 


LEV 5 L 18 


< SEPT 69 ) ;/ 05/360 FORTRAN M ' 

MLER OPTIONS - NAME* MAtN.OPT= 01 .LINECNT»S 8 . SOURCE, EBCDIC. NOL 1 ST .NODECK. LOAD., MAP 

NAME ATEMP" " “ " 


DATE 

i NOED IT# 


72.061/14, 

ID* NO XREF 


> 44*08 


0 

« 


COMP! 


04 ® 
O' 


■s> 

4 ?) 

$ 

* 


■ C 

c 

C LANGUAGE 
C 

c machine 
c 

“C'PURPOSE — 

c 

c 

c 

c 

c 

c ROUTINES REQUIREO 

c 

c 

c 

c 

c 

C TNPUT/DUTPUT 

C 

C 

C 

C 

c 

— C 

c 

C 

c * 

— c "~ 

C ‘ 

c 

c 

c 

c ■ " T 

c 


fortran IV 
IBM 360 

-TO”PREOtCT -TEMPERATURE FOR RAE-D 


ADOYMD 

CROSSV 

NOISPW 

STATIC 


ANGLE 
D!FF 
OROF! X 
01 FFTM 


ARCTAN 
D I NRAO 
PEVERS 
PLUCK 


BASEVR 

DOTTES 

ROTATE 


CMEAN 

MODULO 

PYMDI 


1 ATE " 
i ATE 
1 ATE 
I ATE 
l ATE 
l ATE 
■ 1 ATE - 
1 ATE 
1 ATE 

J AT F 
1 ATE 
1 ATE 
1 ATE 
1 ATE 
l AT F 
1 ATE 
1 ATE 
1 ATF 


WRDC SC 4020 PLOT PACKAGE __ 

FT 04 FC 01 - OISK INPUT < EPH EMERIS ROUTINE! 

FT 05 F 0 CI — PUNCHED CAROS INPUT 

(STATION POS I T I 0 # N * ORBITAL ELEME 
SATELLITE ATTITUDE ^TRANSMITTER 
AND RECEIVER CHARACTERISTICS) 

— FTC 6 F 00 T " PR INTFD OUTPUT" 

(INPUT LISTS* PREDICTED SIGNAL 
STRENGTHS AND NOISES) 

FT 20 F 00 I - MAGNETIC T APE OUTPUT FOR THE 

- - SC 402 C PLOTTER (SIGNAL STRENGTH) 
FT 2 IF 001 - MAGNETIC TAPE INPUT 

- ( TRANSMITTER ANTENNA GAINS) 

PLOTAPE - MAGNETIC TAPE OUTPUT FOR THE CALCOM 1 ATE 
PLOTTER (TEMPERATURE) 


-I" 

2 

3 

4 

5 

6 
7 
B 
9 

10 
1 1 
12 
1 3 
1 4 

15 

16 
17 
1 B 





ISN 

0002'- 

© 

isn 

0 0 3 3 

'3 

ISN 

0 ^ C 4 


TSN 

Z CO 5 " 

•0 

ISN 

0006 

© 

TSN 

C 0 & 7 


ISN 

CC03 


ISN 

0009 

3. 

ISN 

90 IP 


I SN 

con 

9 




— IMPLTCIT—RE'AtTfrB 1 A-"H * 0 -Z ) 


-'DIMENSION LXSATI 3 36) .LYSATC 3361 . LZS A T ( 336 I , LXDOT ( 3 36 ) . 

* LYDOTI 336 ) .LZDOTI 336 ) .SI NMI N! 7 ) , I FSKY < 7 > 

DIMENSION AUNIT13>,EUNIT(3>.GMUNIT(3).1ND(S) 

INTEIMY!' TEST!, TEST 2 . MFSHTR ( 2 66 9 > . MSSHRC . M APFCY .KC . KARY I 336 . 7 > 
INTEGER EPOCMY . EPOCHH. TYPE .ECONS-. 

-S SPACE. CLOCKS. GEOIO. FREQ 

- REAL *4 LXS AT.LY SAT. LZ3AT. LXDOT. LYOOT.LZOOT.LTS ATX. LTSATY.LTS AT Z. 
X LSTAX.LSTAY.LSTAZ.LZHAT.LSLAT.LSLON.LELM^N.LGAIJVA.LEICHT, _ . 

I IWAVE.IRANO.IPANGE. I R SINE .T a 

REAL«A 1 NTGTH, I N TGPH . T AR Y t 3 3 6 . 7 ) ,CTEMP( 7 .71 


REAL*S NAME. IDENT. LIGHT. MEGAHZ 


LOGICAL*l S w l CH • DUMP . R A T I O • T APF 2 1 * T APE 2 2 

LOG ICAL »1 NAME ID! 6 ) .MESH IOC 6 ) . HAP ID I 6 1 . CHECK I 12 ) * CHECKS I 2 ) 


1 ATE 

19 

l ATE 

20 

1 ATF 

21 

) ATE 

22 

1 ATE 

23 

t ATF 

24 

l ATE 

25 

l AT F 

2 6 

1 ATE 

27 

1 ATE 

28 

I ATE 

29 • - 

1 ATE 

30 

1 ATE 

it 

1 ATE 

32 

1 ATE 

33 

! ATF. 

34 

1 ATE 

3 5 ■"* ' 

1 ATE 

36 

1 ATE 

3 7 

1 ATE 

38 

1 ATE 

39 ‘ T " T ~ 

1 ATE 

4 0 

1 ATE 

41 • • 

1 ATE! 

42 

\ ATE 

4 3 - ** * 

1 AT C 

44 

1 ATE 

45 

1 ATE 

4 6 

1 ATE 

47 

l ATE 

40 

~ \ ATE 

49 • 

I ATF 

SO 

1 ATE 

51 

1 ATE 

52 

1 AT E 

53 

1 ATF 

54 

' l AT E 

55 
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ISN CC1 2 

I SN 0613 
ISN 0C14 
ISN 0C l 5 


ISN 0016 ' 


ISN C0J7 


ISN ccia 


ISN 0019 


ISN 0 C 20 
ISN 0C21 


ISN 0 C 22 
ISN OC 23 


ISN 0024 
ISN 0025 
ISN 0 C 26 
ISN 0027 
ISN 0025 
ISN 0029 

ISN CC 30 

ISN CO 3 1 

ISN 0032 
ISN 0033 


LOGICAL*! PUT *PCH . I STP 

EQUIVALENCE < IDENT • N AME I D< 1 > > * < MESHTR< | ) , ME SH 1 D< I I) " ** 

EQUIVALENCE ( TE ST1 .CHECK 1 (!)>♦( TEST 2 .CHECKS ( l ) ) 

EQUIVALENCE <MESHRC(7>, I NO ( 1 ) ) • { MESHRC ( 4 ) , MAPFCV ) • 

X ( MESHRCC 2349) .NOIGROI 

STORAGE SPACE RESERVED FOR ORBIT COMPUTATIONS AND ANTENNA GAINS 

DATA SPACE/33.6/ .nLOCKS/2668/ : ; : 

SEMI MAJOR AXIS OF EARTH .. 

OATA AUNIT /6 378 1 66. , 6370388# • 6 3781 65# / - *** 

FLATTENING- OF EARTH ELLIPSOID 

DATA EUN IT/298#2S.297#0 # 298# 3/ 

SQUARE ROOT OF CM IN UNITS OF EARTH RADIUS CGM/A**3> 

(OR 1/CANONICAL UNITS) 


OATA GMUNIT / #0 0 1 2 394454 71 1 •* 00 1 239 4 1 539254 .# 0 0 1 239443589 / 


DATA DUMP/. FALSE#/. TAPE 21 /.FALSE# /• TAPE 2 2/# F ALSE#/ • CHECK 1 /2H /# 

X CHECK2/2H / - 

DATA MAPID/6HSKYMAP/,NOSKY/1HN/.MEGAHZ/0#DC/ 


SPEEO OF LIGHT IN METERS/SEC 


DATA LIGHT/2#9979250f 8/ 
OATA I I/C/.LL/0/.MI/48/ 


COMMON /INERTL/ ORBELM( 6) .THOOT 
COMMfJN/RF. C6V/ I NTGTM * I N TGPH . ME SHRC ( 18732) 

COMMON /ROTATS/ l TYPE . M A TRX( 21 ) 

COMMON / ST A I O/NAME ( 7 ) . TYPE ( 7 ) • LATO { 7 ) .LATM(?) . LSL AT ( 7 > . LONO( 7 ) * 
S LC1NM( 7) #LSLQN( 7) .LGAINA(7) • LEL M I N, NSTA • LE l GHT ( 7 ) 

COMMON /ORIENT/ LHATSC42 ) ,LZMAT( 3. 7) . ISUN(2 ) . I WAVE# IDAND 


COMMON /MISCEL/EPOCHY.EPnCHM.I NTRVL. I STR T Y# I STRTH# NUMSAT # . 

* L T SA TX.LT5ATY » LTS A TZ .LSTAX (7) *LST AY (7) *L5TAZ(7)» 

* T A . I R^NGE . ! R S I NE 

CO.MMON /NEW/FREQ. PCH ... . 

DATA RD/. 17453292520-1/. BLANK/OH - ~~ / .GEOI 0/0/ * *“ 

NAMELIST /CONTL/ TAPE21 #TAPE22.DUMP 
NAMEL IST/INTGER/ INTGTH . INTGPH 


1 ATE 
1 ATE 
SATE 
l ATE 
I ATE 
l ATE 

•58 

59 

6C 

61 

62 “ 
63 

— A A II 

1 ATE 

O «* 

65 

1 ATE 

66 -* 

1 ATE 

67 

1 ATE 

63 

X a r E 

69 

1 ATE- 

— 70 

1 ATE 

71 

1 ATE 

72 — 

l ATE 

73 

1 AT£ 

74 

1 ATE 

75 

1 ATE 

■76 

1 ATE 

77 

1 A^E 

70 

1 ATE 

79 

I ATE 

89 

1 ATE 

81 

... _ 

1 ATE 
1 ATE 

8i3 

1 ATE 

64 . * 

1 ATE 

85 

1 ATE 

86 

1 ATE 

87 

1 AT F 

** 88 

1 ATE 

89 

1 ATE 

99 

1 ATE 

91 

1 ATE 

9 2 

1 ATE 

93 

1 ATE 

94 

1 ATE 

95 

1 ATE 

96 

1 ATE 

97 

1 ATE 

98 

1 ATE 

99 

1 ATE 

133 " 

1 ATE 

101 

1 ATE 

102 

1 ATE 

103 

1 ATE 

104 " 

1 ATE 

105 

t ATE 

10 6 ~ 

1 ATE 

107 

1 ATE 

108 

1 ATE 

l 0 9 

1 ATE 

110 

1 ATE 

1 1 1 




PACE 093 


ISN 901* 
ISN 9035 


ISN 9036 


ISN 3037 
tSN 0C 39 

tSN 9940 


tSN 3041 
ISN >942 
ISN 0043 
ISN 0944 
I s N ) 9 4 6 
ISN 0047 
ISN 2048 
ISN 0049 


ISN COSO 
ISN 3051 
ISN 90S2 
ISN 9063 


DtflNC CONSTANTS 

THOOT= .72921 15656480-4 
CALC PLOTST ( 1000.41 


-READ CONTROL CARO 


REAOIS . CONTL » 

„„„SFT SKY NO.SE MODEL-. RECEIVER AND TRANSMITTER ANTENNA 
GAIN PATTERN STORAGE TAPE TO LOAD POINT AND 
LOAO TABLES IN MEMORY 

IF< .NOT.TAPE2I » GO TO 5 

READ! 5. INTGER) _ 

REA0(°1^EN0=223.ERR=223) « ME SHTR ( J I . J=1 • BLOCKS » 


1 f p T , y A » SK YM AP^ ‘ THEA) 0 l^T^CONT aTnS^TME^TR ACK I NG RECEIVER 

InTENNA GAINS. RADIO STARS. FLU X 0 E N SITY_ AND_RAD 10 SKY M AP )_ 

DO 1 1-1.6 

CHECK 1 ( 1 l=MAP ID ( I ) 

CHECK2I 1 )=MESH10( I ) 

IE< TEST1.NE.TEST2) GO TO 5 

1 CONTINUE 

00 2 I -l. BLOCKS __ . _ 

2 MESHRC l I I = MES1ITR( I > 

MEGAM*=MAPFCY 

SET LIMIT OF TOTAL NUMBER OF SKY TEMPERATURES TO BE STORED 

latRO/MOIGRO 

- X0UMT=2»I*< Ml) ♦2350 

t SC -KOUNT/DLOCK S 

JFIKnUHT.GT.I 1SCT0L0CKSH ISC=ISCM 


1 ATE 

l 12 

t ATE 

113 

1 ATE 

1 1 A — ' 

l ATE 

115 

l ATE 

l 16 

t ATE 

117 

1 AT E 

1 l 8 

1 ATE 

1 19 

1 ATE 

120 : 

1 ATP 

1 2 l 

1 ATE 

122 

1 ATE 

1 ? 3 

1 ATE 

124 

l ATE 

1 25 

1 ATE 

1 26 

1 ATE 

127 

1 ATE 

12 8 " 

1 ATE 

1 29 

1 ATE 

130 - 

1 ATE 

131 

"1 ATE 

1 32; — 

1 ATE 

1 33 

1 ATE 

134 - 

1 ATE 

1 35 

1 ATE 

1 36 

1 ATE 

1 37 


MAP BLOCK LAYOUT IN UNITS ClF 2 0YTES- 


TABLE 

LOCATION 


NO; OF" ‘ DATA 

CELLS DESCRIPTION 


VARIABLE 

TYPE 


TABLE !0.=*SKYMAP. 

FREQUENCY (MEGAHERTZ) 

RECEIVER GAIN GRID WIDTH 
( .2 PEG. MINIMUM) 

RECEIVER GAIN PRESENCE i NDTC AT0RSI5 ) 
RECEIVER I3FAM WIDTH (DEC. DEG.) 

( 5 CELLS) 

RECEIVER ANTE NNA GA I NS . I OB ) 

(5 x 451 CELLS) ^ 

~ BLANK ’ 

SUN TEMPERATURE (KELVIN DEG.) 
ANGULAR 01AM. OF SUN (DEC. DEG.) 

number of radio Stars (it max.) 

RADIO STARS (10 X 3 CELLS) 

(RIGHT ASC., DECL.. FLUX DENS . ) 


R* A 

> 

m 

163 

P*A 

1 ATE 

164 

1*4 

1 ATE 

165 


1 ATE 

166 

R*4 

1 ATE 

167 
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c 

2 349 

2 

SKY MAP GRID WIDTH (2 DEG* MIN*) 

1*4 

IATE 

i6d 

c 

2 351 

16380 

SKY TEMPERATURES (KELVIN DEG*) 

1*2 

IATE 

169 

c 

tnrjl 

2 

BLANK 

1*2 

— IATE 

170 


V3 


ISN 0955 

ISN 90 56 
ISM Of 5 7 
ISM 00 SB 
ISM C06C ‘ 
I SN 0 C 6 1 

ISM OS S3 

I 5M 0063 
ISN OC 64 


ISM C066 
XSN 3067 


ISM 0060 
tSN 0069 


ISN 0070 


OO 3 1 = 2* ISC 

T 0= ( I- I ) *UL OCK S+ 1 

100=1* BLOCKS .. 

IF{ KOUNT.UT* 100 ) IOO=KOUNT 
•3 P£AD< 2 1 *END=223 *ERR=223) < ME SHRC t J > , Dv! DO>— — 

R E A D( 2 1 *EMD=4 ,ERR=223> (MESHTRI J) . J = l . BLOCKS ) 

GO TO 5 

4 TAPE21=*FALS£* 

5 IF< TAPQ22 ) REWIN022 

— - — RE AO PREDICT ION INTERVAL r 


ISTRTY, ISTRTH ••• START TIME IN YYMMOO HHMMS5 

t STftPY • ISTOPH «•• STOP TIME IN YYMMOO HHMVSS 
INTRVL., «T IMG STEP (IN SECONDS) FOR ORBIT COMPUTATIONS 

RE AO ICOOO. ISTRTY, l STRTH , I STOPY , ISTOPH , T NTRVt 

MM=7*G6400/INTRVL 


-COMPUTE THE DIFFERENCE IN SECONDS BETWEEN START AND STOP 


CALL DIPF( ISTRTY, ISTRTH, I STOPY • I STOPH * I 0 * I SC ) 

“ ITIME= 10*864004- ISC - - 


CALCULATE RIGHT ASCENSION OF GREENWICH AT JAN 0*0 

OF YEAR OF START OF RUM AND ELAPSED DAYS SINCE JANOO 

CALL BASEYR ( I S TRT Y , ! S TR TH , I DO , I SCC .0 AY 1 f 0AY2 f THETOO ) 


C 

c 

c- 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 


-READ STATION CAROS 


NAME *i 
TYPE • . 


I F SKY • 
LATO • 
L ATM • 
LSLAT 
LOMO » 
LOOM . 
LSLON 
LE I GHT 
LELMIN 
LSEMS , 
LGA INA 

MOOIJL . « 
LOSSPW 


• STATION NAME. 

• =1 SATAN SYSTEM 

-2 05-FOOT DISH • * 1 

=3 AT-FOOT OISH 
- 4 16- ELEMENT YAGI ARRAY 
=5 9-ELEMENT VAGI ARRAY 

• SKY NOISE COMPUTATION SUPPRESS CODE 

• DFGPEES OF GEODETIC LATITUDE 

• MINUTESS OF GEODETIC LATITUDE ; 

•• SECONDS OF GEODETIC LATITUDE 

• DEGREES OF FAST LONGITUDE 

• MINUTES OF FAST LONGITUDE 
•• SFCONOS OF EAST LONGITUDE 
•••SPHEROID HEIGHT IN METERS 

• •• MINIMUM ELEVATION ANGLF IN DEGREES * 

•• T R E5HOLD SENSITIVITY (D8M) 

• •• GAIN OF RECEIVING ANTENNA ABOVE ISOTROPIC 

SOURCE (DOW) , » ’ 

•• TELEMETRY RECEIVER DETECTOR ASSYMPTOTE LEVEL (DBM) 
TRANSMISSION LINE POWER LOSS FACTOR 


I ATE 171 
1 AT E 172 
1 A T E 17 3 
1 AT E 174 
1 AT E 175 
-1 ATE 176- 
l A T E 177 
1 ATE ! 7 B 
1 AT E 179 
1 ATE 1 89 ‘ 
1AT.F 101 
l ATE 182" 
1 ATE 183 
I AT E 184- 
IATE 105 
1 AT E 186 
1 ATE 187 
1 ATE 188- 
l AT E 189 
1 AT E 190' 
IATF 191 
1 A T E 192 
l AT E 193 
1AT£ 194" 
IATE 195 
1 ATE 196 

IATE 107 

IATE 198 
IATE 199 

1 ATE 2.90 ' 
IATE 20 l 
IATE 2.0 2 
IATE 293 
I A T F. 204 

IATE 295 
1ATF 206' 
IATE 207 
1 ATE 200 
IATE 209- 
IATE 210 
IATE 211 
IATE 212 
IATE 213 
IATE 214 
IATE 215 
IATE 216 
IATE 217 
'IATE 213- 
1 AT E 219 
IATE 220 
IATE 221 
IATE 222 
IATE 223 
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XTEMP».*RECEIVEP REFERENCE TEMPERATURE (KELVIN DEG*) 

I n X ENT *** TRANSMISSION LINE AMBIENT TEMPERATURE (KELVIN DEG) 
- - ' MODPOL. ••POLARIZATION MODE ( I =L I NEAR , 2=C I RCUL AR > 

DO 1O0 1-1*7 

10 PFAO 1 C 2 r 0 * NAME ( I ) • T YPE < I > *LATO< l ) *L ATM< I > « LSL AT ( l ) *LOND( t ) * 

* LONM ( I ) .LSUONC I > .LGAI NA( I ) ,LE ir,MT( I ) ,UELMIN 


-CHECK' FOR -LAST STA T ! ON ' C A RD - 


ISN 0073 


IF{ NAME< I >.EQ*BLANK) GO TO 200 


-SCREEN OUT INVALID CODE 


IF< TYPEC I )*UT. t ) GO TO 11 

IF ( TYPE ( I 1 *LT® 6) GO TO 20 
U PRINT K41C 
GO TO 10 


-COUNT NUMBER OF STATIONS 


rsN oc3i 


20 NSTA=I 


ISN 0002 
ISN 0l83 


— . COMPUTE SINE SQUARE OF MINtMUM ELEVATION 

OOURNO=LELM IN. 

SIMM 1N( I > = < OS I N< DOUONOORD) } **2 


-PRINT PROGRAM IDENTIFICATION PAGE AND TRACKING STATIONS LIST 


ISN JCBA 
ISN 0035 
ISN 0037 
ISN OC 83 


CALL HEAD IN < I , I STOP Y* I STOPH . I PAGE • KOUNT ) 
IFC XNOCTYPEC 1 ) )*£Q« 1 > GO TO ipO 

PRINT 1CA30* NAMSUI 

K0UNT=KnUNT+2 
ICO CONTtNUE 


PASS BLANK CARO SEPARATING OATA 


ISN 0089 


• READ 1(200* IOENT 


-READ HROUWFR MEAN ELEMENTS AND SATELLITE CONSTANTS 


-IDNO .«« SATELLITE TYPE (l=$PIN TYPE* 2=GRAVtTY STABILIZED) 
l DENT SATELLITE IDENTIFICATION 

EPnCHY,FPOCMH EPOCH IN YYMMDO HHMM5S 

O.. ♦ORBITAL ELEMENT TYPE CODE t BLANK*MEAN*0 = OSCULAT I NG) 

DA ••• SEMI-MAJOR AXIS 
DE ••• ECCENTRICITY 
D I ••• INCLINATION 
DM «.* MEAN ANOMALY 

DP ••• ARGUMENT OF PERIGEE * ; * •* 

DN ••• RIGHT ASCENSION OF ASCENDING NODE 

DO ACCELERATION OF MEAN ANO M ALY < R AD I AN/C ANCNI C AL UNIT**2) 

DT t ME « • o OR AG CONSTANT REFERENCE TIME FROM EPOCH IN SECONDS 
ECONS**. EARTH CONSTANT CODE (0*1.2) 

IPOT* • *ROTATI ON RATE OF SPIN TYPE SATELLITE ( RRVOLUT I ON/M t N. ) 




pace eoa 



C POWER. ..RADIATED POWER (WATTS) IATE 280 

c FCCO ... FREQUENCY IN MEGACYCLES PER SECOND 1 AT E 281 

C BAND* ... BANDWIDTH IN KILOCYCLES PER SECOND . : » ATf ?{>S •— - 

C M INMAX. . .MI N I MUM AND MAXIMUM ANTENNA GAIN (DBW) 1 AT E 283 

C KD ..... DETECTOR CORRECTION DUE TO MODULATION (OBW) 1 ATE 20A 

c 1 AT £ 2(35 

ISN 0090 2C0 NUMSAT=0 1ATE 236 

ISN 0 091 P£ AO 10100. IDENT. EPOCHY .EPOCHH .ECONS.FREO l AT E 28? 

c -• — 1 AT E 289 

c CHECK FOR LAST SATELLITE l*TE 289 

c 1 AT E 290 

ISN 0092 I F( EPOCHY. EO.O ) GO TO S000 . 1ATE 291 

c - 1 ATE 292 

C COUNT SATELLITES 1ATE 293 

C - ‘ — 1 AT E 29A > 

ISN 0094 NUM SAT = NUMSA T+ 1 1ATE 295 

c — — — ' •• - — - ' 1 ATE 296 r- 

C SET EARTH CONSTANT CODE 1ATE 29/ 

ISN 0095 ECONS=ECONS*l IATE 299 — 

C 1 AT E 299 

C SELECT EARTH ELLIPSOID — IATE 30 0 ' 

c l ATE 30 I 

ISN 0596 AE = A UN IT(ECONS) IATE 302 

ISN 0097 F = 1 • DO /EUM IT(ECUNS) IATE 3^3 

ISN 0098 SORTGM=GMUN! T ( ECONS ) 1ATF 304 

C .IATE 305 

C ' — : IATE 3CG 

c LOCATE TRACKING STATIONS ON SELECTEO GEOID IATE 307 

C IATE 30 9 

ISN 0 399 IFIGEOIO.EQ. ECONS) GO TO COO IATE 309 

ISN 0101 GEO 10= ECONS IATE 310 

ISN 0102 E WGL 1 = ( 1 .0-F ) «*2 IATE 311 

ISN 9103 EWGLS0=1 .C-EWGL1 — IATE 312 

ISN 0 1 0 A CALL STATIO ( AE • E WGL 1 *E WGLSO ) IATE 313 

C IATE 314 

c COMPUTE ELAPSED SECONDS BETWEEN EPOCH IATE 315 

C AND START TIME OF RUN IATE 316 

C IATE 317 

ISN 0 105 SOC CALL HKAOINI I STRTY. I STRTH. ISTOPY.I STCPH, INTRVL) IATE 310 

ISN CIC5 CALL DIFF ( EPOCHY .EPOCHH, I STRTY. I STRTH, IOAY. ISEC) IATE 319 

ISN 0117 ISECN= II>4Y*864C0*-I SEC IATE 329 

ISN c ICO PCHs.TRUE. IATE 321 

c 1 ATE 322 

c --SET DATE FOR DROIT CALCULATION'S . 1 AT F 323 

C AND - IATE 3 24 * 

c -—SET JSEC TO ZERO TO KEEP TRACK OF PREDICTION TIME INTERVAL IATE 325 

C 1 ATE 326 

ISN 0 109 ISC=ISCC IATE 32.7 

ISN 0110 JSECNsf. IATE 323 

ISN 0111 FOUNTS t . 1ATC 329 

ISN 0 112 DAYI 1= IDO ' ; IATE 330 

ISN 9113 ISECO=-INTRVL IATE 331 

ISN 0114 INHO=0 IATE 332 

ISN 0115 GO TO 1500 IATE 333 

C IATE 334 

c INCPEMENT TIME BY INTERVAL *ATE 335 
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0 

O 


. c 


-p* 

ro 


ISN C t 16 
ISN 0117 


ISN 0118 


ISN 0120 


ISN 0122 
ISN C 1 23 
ISN 0124 


ISN 0 1 2 1 


ISN C127 


ISN 0128 
ISN 0129 


ISN 013? 


ISN C I 3! 


ISN 0 132 


ISN 0133 


ISN 0 134 
ISN 0135 


ISN 0135 
ISN r 137 
ISN 0138 


c ' • ' - I ATE 336 

1CC0 JSECN-JSECN+tNTRVL 1ATE 337 

ISTP=JSECN.GT.I TIME — — — : 1ATE 3 3 0— 

C * 1 ate 339 

C— CMCCK FOP STOP TIME 1ATE 340 

C l AT £ 341 

IF ( JSECN.GT . T TI ME ) GO TO 2000 - IATE 342 

C . ! ATE 343 

C CHECK FOR STORAGE CAPACITY IATE 344- 

C IATE 345 

IF(KOUNT.GE.SPACE) GO TO 2000 IATE 3^6 

C IATE 347 

C INCREMENT ALL TIMES FOR SATELLITE *• - * ■ IATE 348 *• 

C • l ATE 349 

1250 ISECN= I SECN+tNTRVTT IATE 3 5."? — 

KOUNT«KClUNT+l IATE 351 

I SC* ISC* INTRVL •" * IATE 352 

C IATE 353 

C CHECK FOR NEW DATE IATE 354 

C IATE 355 

IF< ISC.LT. 06400) GO TO l 500 IATE 336“ 

C IATE 357 

C SET NEW OATE " IATE 350 

C IATE 359 

I SC= ISC— 06400 IATE 369 

C IATE 361 

C ; SET ELAPSED DAYS AND OECI MAL * DAYS FOR WIGHT* "AS C OF CRNWCH IATE 362" 

C IATE '363 

DAYH=0AYI 14-1. IATE 364 

I5CC OAY2 = FL0AT( ISO/86400. IATE 365 

C • IATE 366 

C CALCULATE RIGHT ASCENSION OF GREENWICH IATE 367. 

C ' IATE 368 

CALL ROTATE ( THCTCO .DAY I 1 .OA Y2 • THETG ) IATE 369 

C ••••''• IATE 3 7f> 

C CALL OPH I T TO COMPUTE INERTIAL X , Y • 2 . X OOT * YOOT * ZOOT IATE 371 

C l A r E 372 

CALL OR3MSI I SECN • EPOCHY * EPOCHH ) IATE 373 

C ... •— IATE 374 

c TRANSFORM TO EAPTH-FIXEO COORDINATES AND VELOCITIES 1 ATI- 375 

C IATE 376 

CALL OROFI X ( THETG tLXSATl KOUNT) .LYSAT < KOUNT) «L2S AT < KOUNT ) • IATE 377 

X LXDOT(KOUNT) #LYDOT( KOUNT) •LZOOT(KOUNT) ) t A T E 376 

c IATE 379 

GO TO ICOO IATE 380- 

C ... IATE 361 

C 1 ATC 382 

2000 CONTINUE IATE 383 

PCH=. FALSE. ■ IATE 384 

. c IATE 385 

c — — .---COMPUTE SATELLITE VISIBILITY AND TRANSMITTED SIGNAL STRENGTH IATE 306 

C AT EACH STATION IATE 387 

C * 1 ATE 388 

4 2CC PLT=KOUNT.EQ.MM IATE 389 

00 40CC Ksl.NSTA IATE 390 

LL-0 IATE 391 


o 


o 



t SN r 1 39 ■ - 11 = 3 , ••• ■ . IATE 392 

ISN r.|4« ITYPE=tYPE< I > 1 ATE 393 

ISN Cl»l RAT I0=. FALSE. : t ATE 394' 

c 1 AT F 395 

C " CHECK IF RECEIVER ANTENNA PATTERN IS AVAILABLE I ATE 396 

C IATE 397 

ISN 0 142 I F ( IND < I TYPE). ECU I) RATIO*. TRUE. 1 AT E 39ft 

C IATE 399 

C —CHECK IF NOISE PREDICTION IS SUPPRESSED ; IATE 900 

c IATE 401 

ISN C 1 4 4 IF( IFSKYI I ).EQ. NOSKY) R ATI O =» FALSE . IATE 4C2 

C IATE 403 

■ c ........ ...... IATE 404 

C SET PREDICTION DATE ANO CONVERT TIME OF OATE TO SECONDS 1 ATE 405 

c - -• - - — — • IATE 406 

ISN 0146 IYMO=ISTRTY IATE 407 

ISN C 1 4 7 JSEC= I STRTH-40*< I STH TH/ 1 00 ) -24C0 * ( l STRTH /I 0000 > - I NTRVL ' ~ IATE 40S 

ISN 0 14 0 INHniJP=INHO IATE 409 

ISN 0149 I NYEAR= I STRTY/1 OOCO IATE 410 

. C IATE 4 1 1 

C INITIALIZE RIGMTASCENSION AND OECH NAT I ON ' OP SUN IATE 412 

C IATE «13 

ISN 0150 - CALL B A SE YR ( I STRT Y . I STR TH , I 0 . I SC ,0 AY 1 .OA Y2 « TMETGO ) IATE 414 

ISN C151 CALL OPR MSI JSEC .EPOCHY • EPOCHH I IATE 415 

ISN 3152 ISEC=ISECn IATE 416 

C 1 ATE 4 1 7 

■ - C— -« AVF LENGTH CONSTANT ' - IATE 416 

c IATE 419 

ISN 0153 KAVEL=LIGHT/1FRE0*1.0T6) IATE 420 

ISN 0154 IKAVE=AAVEL IATE 421 

c IATE 422 

C IATE 423 

C CHECK EACH COMPUTED SATELLITE TIKE POINT FOR VtSIBILITV IATE 425 

c IATE 426 

ISN 0 1 55 DO 30Cr J=1 .KOUNT IATE JET 

ISN 0 1 56 t 1= I I U IATE 420 

ISN r,|S7 LL=LL»l 1ATF 429 

ISN 0150 JSEC = JSECMNTRVL ' * IATE «30 

ISN 0159 ISEC=ISEC*INTRVL IATE 431 

ISN 0160 1F< JSEC. LT. 064001 GO TO 2100 1A1E «T2 

ISN 0162 CALL ADOYMD t l YMO . 1 ) IATE 433 

ISN 0163 J SEC» JSEC* 8 6400 . IATE 434 

C IATE 435 

c-— CHECK IF ONE HOUR PASSED BY ~ : * ““ IATE 436 

C IATE 437 

ISN 0164 21C-C lESSHR=t t SEC- INHr)UR*360r I /3600 IATE 438 

ISN- 0165 ’ IF ( LESSHR »LE. 0 I GO TO 2400 IATE 439 

C IATE 440 

c COMPUTE R IGHTASCENS ION AND DECLINATION OF SUN IATE 441 

C 1 — IATE 442 

ISN 0167 INHOUR=INHOUR-H IATE 443 

ISN 0160 CALL RIVERS ( J SEC , I HRMSC ) IATE 444 

ISN 0169 CALL B A SEYR < I Y MD . I HRMSC . 10 , I SC . 0 AY l . DAY2 • THETGO > IATE 445 

ISN 0170 CALL ORBMSt ISEC .EPOCHY. EPOCHH) IATE 446 

c IATE 447 


1 SN 0 171 
ISN 0172 
i 6 N 0 171 
ISN 0174 


ISN 0 1 75 


ISN 0176 


— ISN Cl 78 


ISN C 179 
ISN 01 QC 

isn o ie 2 ;* 


ISN C 163 
ISN 0164 
ISN 0165 
ISN 01 66 
ISN 0187 
— ISN 0188 
ISN 0169 

i sn r i9> 
ISN 0191 
ISN 0192 
I SN 0193 
ISN 0194 
ISN p 195 
ISN 0 196 
ISN 0197 


C iate 

C — FORM EARTH-FIXED STATION-SATELLITE VECTOR 1ATE 

2A00 CONTINUE 1ATE 

LT S AT X=L XSATCJ)— LSTAXI I ) - I ATE 

LTSATY=LYSAT< J ) -LSTAY< I) 1ATE 

LTSATZ = LZSAT( J)-LSTA2( I ) IATE 

C 1 ATE 

C— «• — DOT PRODUCT OF STATTOM f $ tOCAlr' VERTICAL— AND IATE 

C STATION-SATELLITE VECTOR IATE 

C 1 ATE 

CALL DOTTED ( L ZHAT < 1 . I ) *LZH A T( 2 , I > *LZHAT ( 3 t X > .LTS ATX .LTSAT V * IATE 

X LTSATZ * IRSINE ) ' IATE 

C IATE 

C —IF IRS I NE "IS NEGATIVE SATELLI TE” I S” BELOW MOP! ZON 1 ATC 

C IAT£ 

IF( IRSINE.LT. 0*0 ) GO TO 2200 ‘ 1 AT ^ 

C IATE 

c COMPUTE ST AT I ON- TO— S A TELL I TE RANGE SQUARED IATE 

c iatf 


c 

c- 

c 


c 

c« 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 


2200 

3CC0 

400C 


PHOSO=LTSATX**2+LTSATY**2*LTSATZ»*2 1 ATE 

IATE 

— CHECK IF SATELLITE IS BELOW MINIMUM ELEVATION IATE 

IATE 

SINSO* IRS! NE**2/RH0S0 ” ’ 1ATF 

IF( SINSO*LT#SINMIN( t )) GO TO 2200 IATE 

IRS lNE = OSORT (SINSa ) IATE 

IATE 

t--1F PASS IS LONGER THAN ONE DAY < SYNCH RON EOUS ) CUT PASS IATE 

*. IATE 

I IF TIME INCREMENT IS LESS THAN ONE MINUTE - CUT PASS WHEN 1ATF 
STORAGE FOR SATELLITE SIGNAL STRENGTH IS FILLED) 1 AT F 

IF( I VIS.CE* 1440 ) GO TO 2200 IATE 

IATE 

— SATELLITE IS VISIBLE* INCREMENT IVIS 1ATF 

1 AT F 

lvts=tvts-fl IATE 

COMPUTE STATION-TO-SATELLITE UNIT VECTOR IATE 

1 A T F 

I R A NGE = D SQR T ( RHOSQ ) IATE 

ltsatx^ltsa tx/irange •• ; iate 

LT0ATY=LTSATY/1 range IATE 

LTSATZr.LTSATZ/IRANGE IATE 

CALL SOUnCEnYMD.JSEC#! *KC) , 1 AT E 

TARY( I I « I ) -TA y ' "■ IATE 

KAR Y(LL « I >=KC 1 ATE 

GO TO 3T 0 0 ~ ' IATE 

TA=0.0 IATE 

T ARY! XI* I )~TA IATE 

CONTINUE ' ' 1 ATE 

CONTINUE IATE 

CALL PRTOUT< TARY#KARV#NSTA«PLT*KOUNT| IATE 

DO 31 1 1- 1 • NST A ‘ • : -• — IATE 

KKsl IATE 


448 

449 

4 50- 

451 

452 

453 

454 

455 

456 

457 

458 

459 

460 

461 

462 

463 

464 

465 
465 
467 
463 
4 69 
4 70 

471 

472 

473 

474 

475 

476 

477 
4 73 
4 79 
480' 
4 fl 1 

482 

483 
434 
4 85 
4 36 
4 37 
4 88 
4Q ; 9 
4 90 
4 9.1 
492 
49.3 
-4 9 4 

49 5 
496 
4 97 
4 93 
499 

50.0 

50.1 
5^2 

50 3 


■ DO 30 I K* 1 * 7 •- - ' 

CTEMP< IK. I 1 >=<>• 

on 30- J l — t * M t ? 3 — — ! 

C TE MP ( 1k» I l)=rAMAXl(T ARY ( KK • t 1 ) »C TE MP (iKtll)) 

KK = K<+1 ' - • ; . 

IF( KK« GT.KOUNT ) GO TO 31 

continue ’ 

L CONTINUE 

WR I TEC 23) C TEMP — 

I F( t STP.OR* oNOT.PLT ) GO TO 5000 

KOUNTsC .1. 

ISECO^ISEC 

’ INM0=INH0UR * ' - 

IF C JSEC.GE. C 86400-INTRVL ) ) CALL ADO Y MO ( I YMD » l ) 

— I STPTY- I YMO — 

GO TO 1250 

l PP I NT 1042^ 

CALL EXIT 

! CALL ENOPLT ....... 

EMOFILE 23 

FORMAT (2(2I6«2X) • I A ) 

FORMAT! A6. 216, 1 X, II . IX, 13) 

FORMAT (A 6* I 1 «2I3»F6«2*1 X»2I3«F6« 2*2X*F5# lit X*F7®2 # IX «F2«0 ) 
FORMAT < 1H0 * 30X* 9 INVALID INPUT CODE*) 

FORMAT ( XH'j » 3C X« • INPUT TAPE TRANSMISSION ERROR*//) 
rORMAT ( 31X * ’RECEIVER ANTENNA GAINS TABLE IS NOT AVAILABLE 
X * FOR STATION *.A6//) 


ISN 0223 




I. IDENTIFICATION 


A. 

Name : 

BLOCK DATA #1 


B. 

Language : 

FORTRAN IV 


C. 

Machine : 

IBM 360 


D. 

Purpose : 

To store alphabetic 
harmonics in common 

constants and zonal 
block . 


II. METHOD 

Not applicable 


HI. CONSTANTS IN COMMON BLOCKS 


Name' 

Variables 

Dimension 

■ Type 

/DAYBLK/ 

MONTH 

(26) 

I 

/ALMN/ 

MONA 

(12) 

I 


MONT 

(12) 

I 

/ANTNA/ 

ANTEN 

(3) 

R*8 

/RECEV/ 

THINTG 

(1) 

R* 4 


PH I NT G 

- (1) 

R*4 


MESHRC 

(1873.2) 

1*2 


Description 

Cummulative days 

Name of each month 

Cummulative days for each 
month 

Name of transmitter antenna 
system. 

Insert constant "zero” to 
these positions. 
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I. IDENTIFICATION 

A. Name: SOURCE 

B. Language: FORTRAN IV 

C. Machine: IBM 360 


D. Purpose: To compute the Antenna-noise temperature 

of a tracking ground station antenna. 

E. Calling Sequence: CALL SOURCE (IYMD , JSEC , I ,KC) 


Name 

Dimension 

Type 

Description 

IYMD 

(1) 

I 

Input-Date (year, month, day) of 
spacecraft's passage over a given 
station. 

JSEC 

(1) 

I 

Input-Time (seconds) of day at 
passage . 

I 

CD 

I 

Input-Tracking station sequence 
number. 

KC . 

CD 

I 

Output-Number indicates the source 
of the temperature. 


F. Common Blocks: 


Name Variable 

Dimension 

Type 

Description 

/ROTATS/ I TYPE 

1 

I 

Input-Select tracking 
system type code.. 

ROTATE 

(3,3) 

R 

. Input-Orientation matrix 
of tracking antenna 

SQUARE 

(l,D 

3 

R 

Input-Vector giving the 
direction of polarization 
in the antenna centered 
coordinated system. 
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Name Variable Dimension Type 

SQUARE 3 R 

( 1 , 2 ) 

SKY * ■ 3 R 

*- • 


Description 

Input-Vector of polarization 
in the earth-fixed topo- 
centric system 

Input -Proj ection of vector 
of polarization in the 
station boresight system. 


,/RECEV/ THINTG 
PHINTG 

BLOCK 

/ORIENT/ WORDS 
ZHAT 

SUNDC 


1 R 

1 R 

18732 1*2 

42 R 


Input-The orientation 
angles increment of 9 
(in deg.) 

Input-The orientation 
angle increment of 4 > 

(in deg.) 

Input-Containing receiver 
antenna gains and sky map. 

Temporary locations 


R Input-Station's vertical 

unit vector 

R Input-Declination of sun 


SUNRA . 1 


R Input-Right ascension of 

sun 


WAVEL 1 


R Input -Wave 1 length of 

transmission 


BANDW 1 


/STAID/ WORDS 2 63 

STAGAN 7 


WORDS 3 2 

/MISCEL/ F0UR1 6 

SATX 1 

SATY 1 

. SATZ 1 


R Input-Effective receiver 

bandwidth 

R Temporary locations 

R Input-Receiver antenna 

gain above isotropic source 

R Temporary locations 

R Temporary locations 

R ) X } Y,Z components of station 

R > moon unit vector in the 

R ) Greenwich-oriented system. 


F0UR2 


21 


R . Temporary locations 
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Name 

Variable 

Dimension 

Type 

Description 


TEM 

1 

R 

Output-Predicted temperature 


F0UR3 

2 

R 

Temporary locations 

/NEW/ 

FREQ 

«§• 

tr 

1 

R 

Input- Frequency (MH^) 


G. Non-System Routines Required: 

REVERS, BASEYR, ROTATE, CROSSV, DOTTED, ARCTAN 


II. METHOD: 

In the computation of total temperature of the moon-to- 
earth telemety link, four contributing temperature sources are 
considered: 

1. Fixed- level source Antenna back lobe temperature 

136 MH Z’ t back = 750 K 
400 MH Z> . T back = 35° K 

2. Variable- level source - 


a) 

Sky 

noises' (galaxy) 

b) 

Sun 

radiation 

c) 

Rad: 

lo stars 


The effects of terrestrial noise sources, the ionosphere 
and troposphere are neglected. 

Front-end thermal noise is created by the internal heat 
supply of electronic components of tracking equipment and it is 
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thus a function of instrumentation system. Receiver nominal 
temperature is therefore a part of tracking station informations, 

• The variable level sources degrad data quality via noise 
temperatures in the antenna main beam and side lobes (Figure 1) 
and their values can be predicted. The method of prediction 
is presented in Section 2. 


III. CONSTANTS AND MESSAGES 

A. Constants: 

Name Value Dimension 

DEGREE 57.2957795 (1) 
RADIAN .01745329 (1) 

TWOPI 6.283185 (1) 

0NE8PI .0397887 (1) 

BOLTZ 1.38 x 10 _23 (1) 

B. Messages: None 


Description 

Angular equivalent of 1 radian. 
Radian equivalent of 1 degree. 


l/8w 

Boltzmann's constant 
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LEVEL 18 ( SEPT 69?;) 


Ln 

to 


C0MOIL FP OPTtONS - NAME 
C 

C NAME 
C 

C PURPOSE 

c 
c 


OS/369 FORTRAN H 
MA!N» OPT «0 1 .LI NECNT=58 • SOURC F « E8C0 ! C « NOL I ST * NODECK 


tOATE 73 . 061 /I 4 . 44.56 


•LOAD.MAP.NOEOIT. I D, NOXREF 
SOURCE 


ISN 

C0C2 

ISM 

c* : v 3 

ISN 

CC }4 

I 3N 

orzs 

ISN 

(. * 96 

ISN 

: v '• 7 

ISN 

S c r: 5 

ISN 

r r S 9 

I SN 

: o is 

I S N 

00 11 

ISN 

•: r, l 2 

ISN 

9 Cl 3 

ISN 

CO 1 4 

ISN 

9C 15 

ISN 

CO 16 

I SN 

00 17 

I SN 

oc i a 

I SN 

R 0 I 0 

ISN 

02 2) 

ISN 

.0 021 

ISN 

0 0 2 2 

ISN 

r ; 23 

ISN 

OS 24 

' ISN 

r?25 

ISN 

0 0 26 

ISN 

ns ?7 

I SN 

3 0 2 8, 

t SN 

G rt 2 9 

I SN 

oo3: 

I SN 

SC 31 

ISN 

CO 32 

ISN 

CC 33 

ISN 

CC 34 

ISN 

n J 35 

ISN 

C C 36 

ISN 

9C 37 

ISN 

CO 35 

I SN 

^0 39 

ISN 

00 4 1 

I SN 

GO 42 


SOURCE 


TO COMPUTE THE TOTAL SYSTEM NOISE POWER OF A TRACKING 
ANTENNA 


SUBROUTINE SOURCE ( I YMD . JSEC • I • KC ) 

DIMENSION ANTFNA< 3) •SIGMTC 3) . STRNOt < 3*tC ) 

I N T EG& R * 2 MESHRC( 451 .5) ,NOIGAL( 1 6380 > .BLOCK, KC 
I NT EGE 9*4 TYPE, FREQ, LIMIT(3) 

PEAL*?* DAY I ,DAY2. THETGC , THETG.PHI 
P C AL * 4 NHAT.TEM.LOriESI 5) 

COMMON /Nt'W/FPPQ .PCM 

COMMON /ROTATS/ I T YPE , R 0 T AT < 3 . 3 ) . SOU AR E ( 3 • 3 ) , SKY < 3 > 

COMMON/RECEV/ THtNTG.PHlNTG, BLCC K( 18732) 

COMMON /OR I E NT/ WORD SI <42).ZH4T(3,7) ,'SUND-C, SUNRA , V AV EL , 8 ANOW 
COMMON /STAlO/WOROS2(63).SrAGAN(7).WORDS3<2) 

COMM.ON/M I SC EL /FOUR 1 ( 6) . SATX . SA T Y, SATZ .FQUR2 (21 ) . TEM . FOU R3 ( 2 ) 
EQUIVALENCE ( BLOCK < 5) .RCGRIO) , C BLOCK ( 17) .LOBES (1 > ) • ( BLOCK ( 27 ), 

X MESHRCf 1 • t ) ) , ( BLOCK! 2283).SUNNOI ),( BLOCK? 22 85) • SUNO I A ) , 

X (BLOCK! 2287). NOSTAR) . ( BLOCK (2289) . STRNO I(i .1)). 

X ( BLOCK ( 2349) .NOIGRO) • ( BLOCK (2351 ) »NOiGAL<l ) ) 

L OG I CAL* 1 PASS. SKIP, RAO, PCH 

EQUIVALENCE ( THSIN.RAStN) ,( ThCQS.RACOS) , (PHSI N, DCS I N) , (PHCOS* 

X DCCOS) , ( SUNEL.STARCL.SUNLON.STOR) 

EQUIVALENCE ( SQUARE ( 1 .1 ) ,ANTENA( 1 ) ) * ( SQUARE C 2 ,2) ,SIGHT( 1 ) ) 


DATA DEGREE /b 7. 29 57 7 95/ , RAD I AN/ . 1 7 4 5 3292 5E- 1 / . TWOP I /6« 283 1 85/ 
DATA PASS/.FALSE./. ONE 8P I / • 0397887/ , BOLT Z / 1 • 38E — 23/ 

DATA L I MIT/1 50, 61 .61 / 

nr at<param)-!0#*mpa«am/ic* > ' ** ’* — 

IDLC U=75-( < l-136>*5/33> 

T E M=r » 0 

RAD=. FALSE. 

KC = 0 

CALL PEVHRS ( J SEC , .1 HRMSC ) 

CALL RA5EYR ( l YMD , I HR MSC . I O • I S • D A Y1 . D AY2 .THETGC > ’ 

CALL ROTATE ( THF TGO , D A Y I ,0A Y2 . THE TG ) 

POT AT( 1. 3 ) = 5 A T K 
POT AT( 2, 3 ) =SATY 
ROTAK 3. 3) =SATZ 

CALL CROSS V ( S A T X , SA T Y . S A TZ . 1 • , 0 o .0 • , ROT AT { 1 . 2 ) . ROT AT < 2 , 2 ) . 

X POTATO, 2)) : 7" ' 

STOR=SOPT(POTATC 1.2) ♦♦2+RaTAT(2,2)**2*R0TATC3.2)*42) 

POTAT( 1.2)=R0TAT( 1 ,2 )/STOR 
POTAT( 2, 2 ) =ROTAT( 2 ,2 ) /STOR 
POTAT( 3. 2 ) =»OTAT(3,2> /STOR 

CALL CROSS V (ROTAT( 1 ,2) ,R0TAT(2.2) .R0TATC3.2) .SATX.SATY.SATZ* 

X ROT AT ( 1,1) ,ROTAT< 2 , i ) .POTATO .1) ) * V 

TEM P=C . 

I F ( PASS ) GO TO 100 

delth=radi an*rcgrid 

DEL PHsR AD I AN4FL OAT ( NO I GRD) 


SOUR 
SOUR 
SOUR 
SOUR 
SOUR 
• SOUR 
— SOUR “ 

SOUR 8 
SOUR 9 
SOUP 19 
SOUR 11 
SOUR 12 
---SOUR "-13“ 
SOUR 
SOUR 
SOUR 
SOUR 
SOUR 
-SOUR- 
SOUR 
SOUR 
SOUP 
SOUR 
SOUR 
-•SOUP 
SOUR 
SOUR 
SOUR 
SOUR 
SOUR 
SOUR 
SOUR 
SOUP 
SOUR 
SOUR 
SOUR 
SOUR 
SOUR 
SOUR 
SOUR 
SOUR 
SOUR 
~ SOUR ' 

SOUP 
SOUR 
SOUR 
SOUR 
SOUR 
‘"SOUR ~ .49 — 
SOUR 50 
SOUR 51 
SOUR 52 
SOUR 53 


1 4 

15 
1 6 

17 

18 

1 9 - 
20 
21 
22 

23 

24 

25 “ 

26 

27 

28 

29 

30 

31 

32 

33 

34 

35 

36 

37 

38 

39 

40 
4 I 
4 2 
43 - 

4 4 

45 

46 
4 7 
4 8 
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ISN 
XSN 
1 SN 
l 5N 
ISN 
I SN 
I ON 

ISN 

ISN 
ISN 
I SN 
r r*,' 

ISN 
I SN 
ISN 

rsN 


CO A3 
00 A A 

CC AO 
OC 50 
T051 
0 0 s 2 
Of 53 
0 0 54. 
or. 55 

vosr> 

0-5 57 
r> ^ 5 8 
0-;59 
COO*' 
8 06 l 


ISN 0062 
r s n r. r 6 1 
ISN CwOA 
ISN 2065 
ISN OC 66 
ISN C'65 
ISN CO 79 
ISN 0071 
(SN 0072 



INCTH=1 

IF( THtNTG*NE»0* > DELTH=R AD I A N* TH I NTG 

IT( THINTCoNF.C* > I NCTHsTMINTC/RCGRI D ... 

IF( PHJNTG.NE.’o. ) DELPH=RAD I AN* PH I NTG 

APE ArOFL TH*OELPH - ' - 

DEL STM=5IN< DEL TH 1 
OFLCTM=COSCDFLTH) 

OCL SPH=S IN( OELPH ) 



LENGT= 1 35/NO IGRD+l 

PASS=*TRUE. ... . 

ICO THSIN=C. 

TGA IN- O. 

SKXPs.THUE* 

2 GAINsHCAT(STAGAN{ I ) ) *DELTH#DELTH 

GA t N-DELTH'C'OELTH 

GO 70 201 • ' 

X01 SKI P=* FALSE# 

200 STOP - THCraS#DELCTM— THSt N*DELSTH — 

IF(K#GT»LIMIT< I TYPE > I GO TO A10 

I F ( (. STC1R fl t E— 6 ) • L T* 0# ) GO TO 410 " — 

THS IN« T»lSIN*OELCTHf THCOS*DELSTH 

THCOS=STOR ■ > 1 ' . . . . ' 

K-K * INCTM 

: GA f N = S T A GAN ( 1 ) ♦MESHRCI Kit TYPE)'”' - 

GAINsMESHPC (K« I TYPE) 

GA t M*nr AT ( GA f N 1 •THS I N*AR£A 
2G1 C IPCLEs-OELPM 
PHS JNa-OEUSPH 
PHCQS= +DELCPH 

3CC C I R CLE = C tPCLC ♦DEL PH r - — - 

IF ( C IPCCK#GE#TWQPI ) GO TO 20C 

STON = PH'»tN*DELCPH+PHCOS*OELSPH ... 

PHCJSt=PMCns*DELCPH-PHSI NAOELSPH 

PHf, iNrSTOP 

ANTENA( 1 ) =THS INAPHCOS 

AN rr.NA ( 2 la TMS!N*PH5f N ' ; - ... 

ANTPNA ( J )=THCOS 
DO 4 00 Mr 1 * 3 

S t G H.T ( M ) s C • 

SK Y{ M ) =c « ' 

DO 4<*C N= I * 3 

ACO SIGHT(M1 = SIGHT(M}+R0TAT(M*NM«ANTENACN> • 1 * 

x C *^ H °^ T ^ e ° T ‘^ AT< 1 ’ 1 > - ZHAT <2'I > -ZHAT<3. [} ,SIGHT(l> ,SIGHT(2), 

! F ( STOP.L r.O. ) GO TO 300 

l THETA=9.’. -ARSINl SIGHT* 3) ) *DEGREE 

CALL A9CTAN ( S I GH T ( 2 ) , S I GHT ( l> , THE T G . PH I | 

tPHI=PHt*OEGPEG ■■■ - . ' - 

!F( tPHI.LT. 1 ) I PH 1 = 1 

M=LENGTv< ( 1PH1-1 t/NOIGROl+I THETA /NO I GRO* 1 
ST09=N0[GAL(M) 

temp=temp+stor*ga:n 

TGA IN=TCA I NfGA I N 


SOUP 

54 

SOUR 

55 

■ SOUR 

56 

SOUR 

57 

SOUR 

50 

SOUR 

59 

SOUR 

60 

SOUR 

61 

SOUR- 

- 6?~ 

SOUR 

63 

SOUR 

64 

SOUP 

65 

SOUR 

66 

SOUP 

67 

-SOUR” 

60" 

SOUR 

69 

SOUR 

70 

SOUR 

71 

SOUR 

72 

SOUR 

73 

SOUR 

74 

SOUR 

75 

SOUR 

7 6 

SOUR 

77 

SOUR 

70 

SOUR 

70 

' SOUR 

O CT - 

SOUR 

01 

SOUR 

02 

SOUR 

83 

SOUR 

84 

SOUR 

85 

sour 

86 - 

SOUR 

87 

SOUR 

88 

SOUR 

8 9 

SOUR 

90 

SOUP 

9 1 

SOUR 

92. 

SOUR * • 

9 3 

SOUR 

94 

SOUR 

95 

SOUR 

96 

SOUR 

97 


90 

99 


SOUP 
SOUP 

sour ion 

SOUR 101 
SOUR 102 
SOUR 103 
SOUR 104— 
SOUR 10'5 
SOUR 106 
SOUR 107 
SOUR 108 
SOUR 109 



PACE C03 




cn 

<C* 


l 



ISN 

0 1C 4 


ISN 

0 106 

— 

ISN 

1137 


tSN 

0 105 


I SN 

C IS 9 


ISN 

C ! 10 


ISN 

5 112 


I SN 

0 113 

... ... 

• ; s n 

5 114 


f SN 

0 115 


I ON 

0 116 


ISN 

C l 1 7 

” 

ISN 

c 1 1 a 


I SN 

Olio 





■' 

* 1 '* ' 

- 

I SN 

3 121 


I S N 

1123 

— 

ISN 

2 124 


ISN 

0 125 



ISN 

0 127 


ISN 

6 1 28 


1 3N 

r i ? 9 


I SN 

C 2 3i* 


ISN* 

0 1 31 


I SN 

1 I 33 


ISN 

C 1 34 


ISN 

0 135 


I SN 

1137 


ISN 

3 1 M 



ISN 

n l 39 


l SN 

Cl 45 


ISN 

3 1A 1 


r sn 

9 142 


ISN 

a i at 


ISN 

0 144 

■ * ■ — . 

1 SN 

0 145 


I SN 

0 1 46 


ISN 

0 147 


ISN 

r- 1 4 5 


ISN 

SIS; 

1 

ISN 

0 151 


I SN 

3 15 3 

' — ■ 

1 SN 

3 1 54 


ISN 

C 1 55 

- -. 

ISN 

0 1 56 


ISN 

C 157 


" ISN 

ft I Sfl 


ISN 

0 159 


ISN 

0 1 60 


ISN 

0 161 


ISN 

0 1 62 


ISN 

0 2 6 4 


IFISKIP) GO TO 101 
GO TO 300 

A 1 0 T fMPt: T (7MP/ TG A ! N - — — ; — 

KC-KC+1 # 

STEMPnr.g 

IF(.3UNSiOI*EO.C#*l GO TO 500 

Sunlon=sunma-thetg - 

OC S I N - 5 I N ( SUN DC ) 

DC CO 5= COS ( SUN DC ) — 

PAS S’ja S INC SUNLON) 

RACOf»=COS( 3UNL.0N) 

ANTENAC l)snCCOS*PACOS 
ANTFNA ( 2 ) =nCCnr.*OASI N , 

ANTETNA ( 3 >=OC3tN 

CALL CGTTSDtZHATC 1 . I ) *ZHATC2.I ) • ZHATOtM t ANTENAC I f.ANTENA(2 V. 

S ANTENAC 3 > •SlfNcL) 

IFC SU.Nfll. .LT.C. > GO TO 500 

CALL dotted C SAT K*5ATY*SATZ* ANTENAC n * AN TE N A < 2 > . ANT ENAC 3 I » SUNEL ) 

SUNEL = ARCOS ( SUN EL ) *0 EG WEE 
IFC SUNEL.GT* 45. 0 I GO TO 500 

' * KC= KC ♦ 2 

IKVSsSUNr.L/’RCGRtD+l 
MnSHS=*.‘.-3HPC < INVS, I TYPE) 

SGA iNshifAT ( MESHS > 

FOR NARROW UEAV5 USE 1/2 POWER OEAMWIDTH 

I F{ LOHFSC I TYPE ) .LT. 15*) TGA I Na< LOBE SC I TYPE > ♦«ADt AN) **2 


S TE Mf*s? ( SGA I N* SUNNO I * CSUNDI A4RADIAN)$*2)/TGAIN 
SCO PTFMsn # 0 

JFC NOSTAR*LE*G 1 GO TO 610 
no 60.1 McIiNUSTAP 
ST0P = 5TON0l C t .M) ♦PAD l AN- THETG 

PAS IN3 5|H( STOP ) - “ ■ 

AAC07.S cost r.rou ) 

STOOsSTHNOl C 2 « M ) “RADIAN 
OCSlNsSJ N< STOP 2 
DCC03=C0 >( STORI 
A N T f A C l ) = u c c n 5 A r. o s 
ant rr.NA ( .» ) :nccn*;*RA3i n 

ANTENAC 3)=DC5IN 

CAUL Q n T r [£ 0 (2HAT( l . I >> ZHAT ( 2 . ! ) , ZHAT( 3. I ) . ANTENA(t) ,ANTENA( 2 », 

X ANTENAC 0> tSTARELI 
1 PC STA»>; L*l.T,C« ) GO TO 600 

CALL DOTTED (SATX*SATY,5ATZ#ANTENACl >* ANTENAC 2 I • ANTENAC 3 J . ST AREL > 

TF( STARL t L.LT.C.9) GO TO 6C0 ; - 

PAD‘.T PUT® 

I THE TA.= ARCOS C r.TAPEL ) ^DEGREE .. .. 

K= T THE! A/RCGRIO +1 

GA r N = Mf- SHRC (K ♦ I TYPE ) 

P TEMPS HEAT ( GA (N ) *{ WAV EL *<’2 ) 4 ST R NO I C 3 . M) «*< ONEQPl/BOLTZ ) 

600 PTEM°al)*0 

ec i ptem=ptc ! urtfmp .... 

t2C CONTINUE; 

IF(RAD) KC-KC+4 

610 TEM=T£«PfSTEMPtRTEH+ IBL(FPEO) 


SOUR 

SOUR 

• SOUR 
SOUR 
SOUR 
SOUR 
SOUR 
SOUR 

* SOUR 
SOUP 

SOUR 

SOUR 

SOUR 

SOUR 

"SOUR- 

SOUR 

SOUR 

SOUR 

SOUR 

SOUR 

"SOUR 

SOUR 

SOUR 

SOUR 

SOUR 

SOUR 

SOUR 

SOUR' 

SOUR 

SOUR 

SOUR 

SOUR 

SOUR 

SOUP 

SOUR 

SOUR 

SOUP 

SOUR 

5nuii 

SOUR 

SOUR 

SOUP 

SOUR 

SOUR 

SOUR 

SOUR 

. SOUR 
SOUP 
SOUR 
SOUR 
SOUR • 
SOUR 
SOUR 
SOUR 
SOUP 
SOUR 


1 10 

in 

t 1 2- 
1 X 3 
1 14 
l l 5 
1 1 6 
1 1 7 
"110- 
1 l 9 
1 20 - 
121 
I 22 ” 
1 23 
1 24*“ 
1 25 
126 
127 
123 
1 29 
1 30“ 
1 31 
1 3 2 

133 

134 
1 35 

1 36- 
1 37 
1 30 
1 30 
140 

1 4 l 

142 - 

143 
1 44 
145 
1 4 6 

1 4 7 
1 4Q 
1 40 
153 
151 
l 52 
l 5 3 
1 54 “ 
1 55 
1 S 6 
1 57 
1 58 
1 59 
1 60 *- 
161 
162 
1 63 

164 

165 





n 

O 

*} 

o 

o 

O 

O 

O 

o 

o 

o 


ISN 3165 
l SN 0166 


RETURN 

FND 


SOUR 166 
SOUR. 167 


3.4 DESCRIPTIONS OF NEW ROUTINES 

I. IDENTIFICATION 


A. Name: EPHEM 

B. Language: FORTRAN IV 

C. Machine: IBM 360 

D. Purpose: 1) Read Lunar 5 Solar Ephemerides and 

Nutation in Right Ascension and the 
Lunar Libration Matrix. 


2) Interpolate the Data Using Fifth 
Order Everett Scheme. 


E . Calling Sequence: CALL EPHEM (TSEC, AO, IYMDE, 

I HME , SECE) 


Name 

Dimension 

Type 


Description 

TSEC 

CD 

R*8 

Input-Time at 
in days from 
year for the 

; which data is desired 
Jan 0.. 0 of the reference 
ARC in seconds. 

AO 

(24) 

R*8 

Output-1) If 
AO 

’ONLYEQ' is true; 
(1) = Nutation. 


2) If 'ONLYEQ ' is false; 

AO (1) -AO [3) - Unit vector 
to moon in 
true of date 
coordinates. 

AO (4) - Range to moon in 
meters . 

AO(5)-AO(7) - Unit vector 
to sun in 
true of date 
coordinates . 

AO (8) - Range to sun in 
meters. 

AO(9) - Nutation. 

IYMDE (1) 1*4 Input -Year, month and day of Epoch 

time 


• 55 



Name 

Dimension 

T yP e 

Description 

IMME 

CD 

1*4 

Input-Hour and minutes of Epoch 
time 

SECE 

CD 

1*4 

Input-Seconds of Epoch time 


F. Common Blocks: None 


G. Non-System Routines Required: DIFFTM 

H. References: • NONAME SYSTEMS DESCRIPTION* 

Vol. I Sec. 2.3.5, 2.4 
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COMPILER 
C * 


ISN 

ISN 

15N 
tSN 
ISN 
I SN 

r on 

I SN 
ISN 
ISN 
I SN 




UP 




Ln <£, 


c 0 

. 


coc-2 

3DM 

00-)4 

3 O'; 5 

?CC6 

0 Z~7 

o * j r o 
0 o 0 9 
CC 1 ■) 
nr, n 
0 0 12 


tSN 
I SN 
ION 

I “iN 

1 SN 
I SN 
ISN 
ION 
I ON 
I ON 
ION 
I ON 
ION 
ISN 
I ON 
I SN 
ION 
ION 
! ON 
I ON 
ISN 

ISN 

r on 
ISN 
I SN 

ION 
ISN 
ISN 
ISN 
I ON 
I SN 
ISN 
I SN 
I SN 
ISN 
ISN. 
ISN 
I SN 
ISN 


CO! 3 
?C 1 4 

oo ; s 
> o i o 
o o 1 ‘1 

: o 19 
OC 2 1 
C 0 22 
?C 23 

.-i 

0 2 26 
0.0 ? 7 
0 020 
CO 29 
r o 

0 P .1 1 
90 32 
■:c 3 3 
CO 34 
O 5 3*3 
?C37 
•0 9 3 8 
)Z 39 
C 3 43 
C 4 l 
«; 9 4 2 
CO 4 3 
0 0 44 
0 045 
*0 4 6 
0 042 
CC 49 
or 50 
3052 
C953 
0CS4 
CO 55 
9 0 06 
0057 


OPTIONS - NAMES MAIN,nPT=01 • L f NEC NT =50 . SOURCE ♦ E0CD I C • NOL I ST • NODECK • LOAD . 

* * EPHEM 

SU3POUTINE EPH£M(TSEC,AO .IYMDE. IHME.SECE) FPHF~ 

PEAL*0 AO! 24) .DAY, PLANET! 3.3. 3.5) .DAYR.FO .F2.F4.FS0 . EPHE 

• FS2.FS4,F2S*F4S,FCS f TSE-C ,OA Y l . MOON ( 3 . 3 . 1 7 ) , F ACTOR .S . TYMOLD EPHE 

P F AL * 0 nUFMl (51 ) ,RUFM2( 5 1 ) • B UF M 3 ( 5 1 ) .BUFPI C 27 > ,0UFP2<54 ) • BUFP3 C 54 ) EPHF 
MF.AL40 OUFM( 1S3) •BUFP(‘135> EPHE 

PFAL NUT (51) • EPHE 

" LOS IC AL NOT I ST/.FALSE./ “ EPHE “ 

INTEGER DA TP/4/ EPHF 

EQU IVALENCE! BUFM( 1 ) • MOON ( 1.1*1)) . (BUFP( 1 ) .PLANET ( 1 * l • 1 • 1 ) > EPHE 

EQUIVALENCE (FSO.S) ephe 

EQUtVALENCE<NUT( 1 ) .BUFMl ( 1 ) ) ' " 

. EOUl VALENCE! 0UFM1 { l > ,0UFM( 1 ) ) • C SUFM2< I > .BUFMi52) ) • 

' • CBUFM3C I ) . f)UFM( 103) ) , OUFPl ( I ) • 0UFP f |*> 1 * ( BUFP2 ( 1 > .BUFPI 28 ) )• — 

• ( BUFP3C l ) .RUFPC 82) ) 

DATA FACrON/U215C3730l 64520-02/ 

DATA MFO/5/. TYMOLD/l • D50/.DAYR/9999. 00/ 

T2( S)=< S*»2-l .DC) *S/6.DO 
I K< NOT 1ST) SO TO 6 
NO r l!iTs, TRUE. 

tF( DAYft.GE.9999.OD0 ) GO TO 5 

CALL DIFFTMC IYMO. IHM.SEC . I YMDE • IHME , SECE.TMlN) 

DAYR=TMIN/ I . 4403 
GO TO 5 

i IF ( TSFC.FO. TYMOLD) RETURN 

5 OAY=T5EC/S. 6404 

IF( CAY-0AV.9 ) 10, 70 .59 
PE WIND DA TP 

READ! DATP.END-*20C ) I YMO . I HM • SEC • NUT * QUFP1 
PEA O(T) A rP.FMO--2O0 ) fJUFMl 
PEA0!DATP.END=2Q0IBUFm2 

READ! DA TP ,END=200 ) l JUFM3 — ~ 

READ! DA TP, CND«2CC )nurp2 
Rfc'AO( n A TP. l:.NO = 20 9 MHJFP3 
|F( IYMD.CO.O) GO TO 209 

CALL DIFFTM! IYMD* I HM , SEC * | YMDE • I HME . SECE t TM I N > 

OAYRvTM in/ l .44 0 3 

IF( OAY-DAYP ) 3 OC • 7C .50 — 

READ! DA TP ,CND=20C ) t YMO. I MM. SEC • NUT • BUFPI 
PF ADC DA r P, LN0=209 > BUF Ml 
PE A D< DA TP • IT Nf>-2 r 9 )BUFM2 
READ! OATP,rNI>=2CC M3UFM3 
PE AD! DA TP, END- 2')0 ) BUFP2 

READ! DATP,F NDs.?C9 )HUrP3 " 

IF! [YMD.FO.O » GC1 TO 200 

CALL OIFFTM! t YMD. I HM . S.EC . I YMDE . I HME . SECE^TM I N ) 

DATO = TM IN/ 1 « 4 4D 3 

IF{ DAY. GT.DAYRf 8.00 ) GO TO 30 

INDCX= IDINT! (DA Y-OAYR ) *2.00 ) 

DAY l = DA YR • 5DO*DFLOAT( INDEX) 

INDEX= INDEX+1 
S=( 0 AY- DAY 1 ) *2*000 
FC S~ 1 . DO — S 
'F2S = F2! FOS) 


MAP.NOFOIT. ID.nDxREF 


2 

3 

4 

5 

6 

— 7 

8 

9 

10 


1C 


20 

3C 


50 

70 


EPHE 

1 1 

EPHE 

12 

EPHE — 

1 3 

EPHE 

1 4 

EPHE 

15 

EPHE 

16 

EPHE 

17 

EPHE 

I 0 

EPHF — 

19— 

EPHE 

20 

EPHE 

21 

FPHF 

22 

EPHE 

23 

EPHE 

24 

EPHE " 

25 

EPHE 

26 

EPHE 

27 

EPHE 

28 

Ephe 

29 

EPHE 

30 

EPHE 

3! 

EPHE 

32 

EPHF 

33 

EPHE 

34 

EPHE 

35 

EPHE 

36 

EPHF 

37 

EPHE 

3Q 

EPHE 

39 

EPHE. 

4 C 

EPHF 

4 1 

EPHF 

42 

FPHF 

43 

FPHE 

44 

EPHF 

45 

EPHE 

46 

FPHF 

47 

EPHF 

40 

EPHE' 

49 

EPHE 

50 

EPHE 

51 ** - 

EPHE 

52 

EPHE 

53 




■PACE 002 


F4 S=<FES**2-4.00 ) 4F2S0.05DO - 

FS2-F2I S> 

FS4 = < S*. 2-4, O0)*FS2*, 0500 1 

DO 90 I I = 1 , 3 

A?< I I ) ^MOONI 1 , I I , I NOEX) *F0S4 MOON! 2, | t .INDEX ) *F2S.M00N< 3, 1 1 , 
•Fas. MOON! 1 . 1 1 ,'INDE X. 1 ) .FSO.MOONt 2 . I t . INDEX* I ) *FS2. 

.MOrjN{ 3, I I , INDEX+n»FS4 ■ . 

! CONTINUE 

INDEX= 10 t NT ( < OAY-DAYR )*. 25>— - 

DAY 1 =PAYR.4,DC *DFLOAT< INDEX) 

IND£X= INDEX+1 - 

■ S=( OAY-DAYl )*.25D0 

Ff.S=l ,09-S • ' 

F2S=F2(F0S) 

F4 

F52=F2( 5) 

FS4 = < S.3-4 . 00 > *F S2* • 0500 

DO 95 J-l.NEQ 

I I = 3 . J ' . 

00 95 [ = 1 .3 

1 I I = I I 4 I • • . V. 

Af( t I I ) =PLANET{ 1,1,1 NDEX.Jl.FOS. PL A NE T ( 2 * I • INDEX , J) *F2S. 
.PLANET ( 3 i I . INDEX, J) .F4S.PL ANET( 1 , f , 1 NDEX.l . J ) *FSO . 

.PLANET (2. I . INDEX. 1, J I.FS 2. PLANE T ( 3 , 1 . I NDEX.l. J)*FS4 
CONTINUE 
DO 96 J - 1 , 3 

ACt J.3 )SA0( J»3>«PAC'T0R*A0( J) 

DO 97 J=2,NEQ 
I=3*J 

DO 97 L - I , 3 

LL=I.L 

AC( LL I =A0 ILL )+A0(L*3) 

tymold.tsec - — 

RETURN 

C WR t TE< 6, 10CC ) I YMD.IHM.SEC . .. 

STOP 12345 

3 WP I TE( 6, 23CO II YMO, IHM.SEC 
STOP 54321 

r FC7RMAT( ' 1PP0GRAM TERM I N ATEO '/'Ol NSUFF I Cl ENT EPHEMERIS DATA*, 
1 ‘CLAST DATA POINT', 18,16. Ffl. 5/) 

: FORMATt • 1PR0GRAH TER M I NA TED • / < 0 I NSUFF I Cl ENT EPHEMERIS DATA*/ 
l • ‘OFIRST DATA PO I NT * , I 8 , I 6 ,F0 , 5/ ) 

END 


- EPHF- 
ephe 

“56 

57 

*EPH£ 

58 

EPHE 

59 

EPHE 

60 - 

EPHE 

61 

ephe 

62 

EPHE 

63 

EPHE 

- 64 - 

EPHE 

65 

EPHE 

66 — 

EPHE 

67 

EPHE 

69 

EPHE 

70 ' * 

EPHE 

71 

EPHE 

72 “ 

EPHE 

73 

EPHE " 

- 74r 

EPHE 

75 

EPHE 

76 

EPHE 

77 

EPHE 

78 

EPHE 

79 

epmf 

80 

EPHE 

81 

EPHE 

82 

EPHE 

8 3 

EPHE 

84 

EPHE 

85 

EPHE ‘ 

86 

EPHE 

87 

EPHE 

88 

EPHE 

09 

EPHE 

90 

EPHE 

91 

EPHE 

92 

EPHE 

93 

EPHE 

94 

EPHE 

95 

EPHE 

96 






I. IDENTIFICATION 

A. Name: HEADIN' 

B. Language: FORTRAN IV 

C. Machine: IBM 360 

D. Purpose: To print the program identification page, 

tracking stations list and radio stars list. 



E. Calling Sequence 

: CALL 

HEADIN (ISTRTY, ISTRTH, 

ISTOPY, ISTOPH - , INTRVL) 

Name 

Dimension 

T y£. e . 


Description " r : 

ISTRTY 

(1) 

1*4 

Input- 
of the 

Date (year, month, day) 
predicted start. . A 

ISTRTH 

(1) 

1*4 

Input- 
of the 

Date (hour , minute , second) 
predicted start. 

ISTOPY 

(1) 

1*4 

Input - 
of the 

Date (year, month, day) 
predicted stop. 

ISTOPH 

(1) 

1*4 

Input- 
of the 

Date (hour, minute, second) 
predicted stop. 

INTRVL 

(D 

1*4 

Input - 

Tracking interval (seconds) 


F. Common' Blocks : 



Name 

Variable 

Dimens ion 

Type 

Description 

/ALMN/ 

MONA 

(12) 

1*4 

Input-Name of each month 


MONT 

(12) 

1*4 

Input- Cummulative days of 
one year. 

/STAID/ 

see 

main routine ATEMP 


/NEW/ 

FREQ 

(1) 

1*4 

Input-Frequency of satellite 
transmitter (MH-) 


59 



Name Variab le Dimension Type Description 

/ANTNA/ ANTEN (3) R*8 Input-Tracking system type 

name. 


G. Non-System Routines Required: RYMDI 


6.0 


trvfL 10 I 


»rPT oo I 


OS/369 FCJTITRAN H 


DATE 72 . 061 / 14 * 44.40 


<!>' 


*4 


& 


if 






JO 

£ 

a 


COMPILER 

C 


OPT ION’S 


NAME = 


!$N 
i :.s 

! 3 N 

I 5N 
I SN 
I SN 


0C92 
? 0 0 3 
7 On a, 
vG 5 5 
4 G 0 6 
O'. 0 7 


i j • < 

r sn 
i r. n 
I SN 


:>7 


ION 

ION 

ISN 


C-C I 2 
OC 1 3 
00 i 4 


ISN 

I $N 

IS N 
ION 
I ' N 
ISN 
I ON 
ISN 


0015 
JO 1 6 

or i 7 
30 10 
0 r i 9 
OC 20 
3-21 
0622 


ION 


00 2 3 
0 0 24 

cr 25 


I 0 N 
i sn 
:$n 
ion 
I ON 
ISN 
I ON 


CO 26 

T 0 

3:20 

0C20 


0 0.3 I 
0 0 32 


ISN 0C33 


ISN 

rsN 


MAlN,CPTrCl ,LlNECNT-=5B.SOUFCC.EOCDXC.NOLXST*NODeCK.LOAO .MAP. NO ED IT, IO.NOXRcF 

HEAOIN 

SUBROUTINE KfeAOINt ISTRTY • 1 S 7RTM . ESTOPVi! STOP*. INTRVLl 

O I ENS I ON FLUXE S ( 30 ) ..DM tf D f 3 ) 

INTECEP*2 MESH 
INTEGER** TYPE. FREQ 
P E AL * 4 INTCiTH* INTC^M 
PE AL*0 NAME, ANTf;N* PONAMC (5) 

LOGICAL*! PCM . — 


5 0 35 


ISN 

ISN 


5 C* 3 5 
:037 


COMMON 
comm :n 
COMMON 
S LONM 

COMMON 
’ COMMON 

EC’J i valence: ( 

5 SUMO I A > • ( ME 
5 ( M*:SH( 2-7 AO >• 


/ A L M N/MDM A { 1 2 1 , MONT( 121 
/fil-CLV/lNTGIM, INTCPH.MESH( 16732) 

/ S r A I 0 / N A M E < 7) * TYPE ( 7>,LATDI7> .LATMt 7 ) ♦ SLAT ( 7 > » LOND ( 7 > * 
7) , SL ON ( 7) .GAIN' 75 . ELM I IM.NSTA .EIGHT ( 7 ) 

/fiFW/."RPO.PCH 

/ANTN4/ANTEN t 3 ) 

SH( 17 ) , QMWOC 1 ) I<* ( MESH (2203) • TEMP ) • ( MESH( 2285 > • 

( 22:37 ) .NUMERO) . ( MESH (2289 > .FLUXES! X > ) * 

V 105) 


DATA 
• CALL 
CALL 
CALL 
CALL 
m :ht= 


R ON A ME/ • CASS I OPEC Y GNUS ATAURUS ACENTA’JRUV I RGO A*/ 
RY.MO I ( I STR T Y , [ Y 1 . I M l , I 0 l ) 

RYMOl ( I STR TH * I Ml , INI, I SI ) 

RYMOK I STOP Y, IY2.IM2.I02) 

RYMD I (I STQPH. I HZ .*!N2. I S2) 

:TRVL/60 


2C 

99 

ICC 


WRITE < 6 * 99) ; 

WPI 71: ( 6 . in.: ) MQNA ( 1 M 1 ) . I D i . X Y l , I H 1 . INI .MONA < I M2 > . I D 2 ♦ I Y 2 . I H2 . 

• WRI TT( Cm If l ) 

DO 1C I*1.N5TA 

1C Wt? [ TEC 6 . 102 ) NAME ( 1 ) . ANTEN( TYPE < I ) ) * L ATO ( I 1 .L AT M ( I ) .SLAT ( X ) » 

Z L O'lO ( I ) ,LONM( I ) , SLONC I ) ,UMWO( TYPE ( I 5) » G A If « { X ) 

XR I If < 6 ♦ t 0 4 ) Tfc MP , SUMO I A 

on 2 ^ K 1 , NU-XE 90 

j=( <- n * if i 

WM I TE < 6 # 1 s; 3 ) PONAMF. ( K ) . FLUXE S( J > .FLUXES< J+ 1 ) .FLUXES ( J*2 ) 

FORMAT ( Jill .////////// ) 

FO*)MAT{ 35Xf MJAE-B TRACKING ANTENNA TEMPERATURE' // 52 X « 

•PER 100 CE COVERAGE * « //» 30X . A4 , \ X . 12 , • . • , • I 9 * , I 2 . I X i 

12** : ' » I 2 » 3 X. * 10* . 3X.A4 , 1. X . X 2 * * • ' • » 1 9 * • I 2 . 1 X . I 2' # • X •» 

////4$x. « AT* , IX. 14, lx. 'MINUTES* INTERVAL *.////.AOX, 

• FPCOUf.NCY = ♦ . 1 X. I 3.1 X. •MHZ * ) 

FORMA T ( 1M 1 , 49 X • * TRACKING STATIONS * ./// . 2X . • ST AT ION* , 4'X • 

S •ANTENNA 5X. «LA T I TtJOF *, 1 C X i- 'LONGITUDE », 10X . 1 ANTENNA •• 

Z- 1GX«* ANTENNA • .//3X, * I DENT. • , 

Z 5 x , * T YRE • . 

• S7X. 'OEG-MIN-Sc.C • t7X« ’OHG-MI N-SFC • ,3X.«MAIN LOOE BEAM tfIOTH«* 

Z 12X , «GA !N* , //*64< , 'OEG. • .23 X , » ) 

1 C 2 rOHMAT ( /y , 2X . AO . A6.5X. I 3 .2X , I 2 . 2X ,F5,Z .5X . 3 3 • 1 X. t2 r 1 X. 

S F5«2. 1 1X.F5.2. 18X.F5.2 .1 I X.F5*2) 

121 FORM AT < ////S 6 X t ’PAD I 0 S T AR S * , / // . 3 3 x . • ST AP NAME* »7X ' “ 

5 'RIGHT ASCENSION* . 6 X. 'DECLINATION* .l.X.'FLUX DENSITY** 

* //.SAX. *DEG* •* 16X • *OEG. *« 1 OX. ‘WATTS. .ME TE RS-2 * HZ- 1 • ) 

1 G 3 FORMAT { //.3CX.A12. X 4 X . F 7 « 3 . 1 3 X • F 0 • 2 1 1 4 X • E9 » 2 ) 

1 C 4 FORMAT! ini'. //« I OX. 5 TEMPERA TUBE FOR QUIET SUN IDEAL MODEL « ff *E9o2 
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1C 1 


HEAD 

■“ t 
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2 
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4 

HEAD 

5 

HEAD 
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7 

HEAD 
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IS 
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27 

| 4 C A D 
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37 
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HEAD 
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I. IDENTIFICATION 


A. 

Name : 

OUTPUT 

B. 

Language 

4 

lb. 

: FORTRAN IV 

C. 

Machine : 

IBM 360 

D. 

Purpose : 

To print out the predicted antenna-noise 
temperature results. 

E. 

Calling 

Sequence: CALL OUTPUT (TSEC, NAME, TA, 

KC, TIME) 


Name 

Dimension 

Type 

Description 

TSEC 

(1) 

1*4 

Input-Total seconds from starting 
time. 

NAME 

CD 

R*8 

Input-Name of tracking station. 

TA 

CD 

R*4 

Input-Predicted temperature. 

KC 

CD • 

1*4 

Input-Number indicate the temperature 
sources. 

TIME 

CD 

1*1 

Input- "TRUE". , print out the predicted 
time 

"FALSE", the predicted time 
is not printed out. 


F. Common Block: None 

G. Non-System Routines Required: None 
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LEVEL 18 l SEPT 69 1 


OS/369 rOPTRAN W 

WMN*t)PT=i)l *t_INECNT*58 # $0URCH*E8CDIC *NOLTST # 


OATH 72*061/14.44*48 


C* 

CrJ 


ISN 
t *’.M 
t ‘i 
1 ■ . ■ ( 

l 

7SN 
1 SN 
ISM 
: sn • 
ISN 
* ISM 
ISN 
ISN 
1 SN 
ISN 
ISN 

ISN 

I VI 
ISN 
ISN 
ISN 
ISN 
ISN 
ISN 
ISN 
ISN 
ISN 
ISN 
ISN 
ISN 
ISN 
ISN 
I SN 
ISN 


0*02 

>*: c. 5 

0 *' * 


1 0 0 f » 

f.*:cr 

orr*8 

0070 

oo n 
•: in 
0 012 

0 0 3 4 
0*15 ’ 
V 0 1 O 

001 r 

CO 10 
'} 3 29 

0 0 22 
0i21 
0 0 25 
0" 27 
0 0 25 
0 0 29 
CO ?0 
0031 
0 0 3? 
0 0 3 3 
00 34 
00 35 
C036 
C 0 3 7 
CC. 38 
0 0 39 


COMPILER OPTIONS - "NAME* 

C . 

SUOPOUT1NE OUTPUT < T SEC *NAM£ * TA **0* T 1 ME ) "* 

I rJ T £ GF P * 2 «. 

I HU. < ,1 V + 4 11 11 . hi >« ♦ r.S * T Sr. C . TOL 0 / - 1 ✓ 

U I «:l . * ‘-t At. 1 *At..£ • Al.3* ill.* 

l_Or.lCAt.4S TIME 


30 


XI 


12 


21 


10 


data At. 1 /'OH SKY T£M»/f AL2/8HSUN T£M*/ f AL3/8H RAOtO* 

‘ OAT A m_K/* * A 

IND=1 

DO 30 1 = 1*3 ’ * 

0(11 =B LK 

I El KC*LT* 41 GO TO 11 ~~ . — J 

KC=KC-4 . 

Ct IN0)=AL3 •• : 

1T4D=XNDM 

TF(KC*LT*2) 00. TO 12 ' ' ’ 

XC=KC-2 

C< I NO) =AL2 ’ *' \ 

T NO-’ l NO 4 1 

IF! TIME-) GO TO 21 
I Ft T SFC • EQ * TOLD ) GO TO 10 
HH=TSFC/36 r 0. 

VM= ( TSEC-HH* 3600/60 
SS= Tsrc-»HH*36P0-*4M460 

Vf? I TE I 6*100) HH#MM,SS*TA*NAME#T> ■ 

TOi. 0 = T SEC 
T IMF. » .FALSE# 

return * 

WRITE! 6* 1021 
TOLOsTSEC 
RETURN 


TA * NAME .G 


LOAD* MAP# NOE D IT# X D* NO XREF 
OUTPUT 



DUTP" 

OUTP 

OUTP 

ti«M P 

ourp 

OUTP 

2 troiR 

r 


6 


OUTP - 

~ T 



OUTP 

a 



OUTP 

9 - - *• 


■ 

OUTP 

10 



OUTP 

1 1 



OUTP 

1 2 



OUTP 

- 1 3 



OUTP 

14 



OUTP 

15 



OUTP 

16 



OUTP 

17 - - 


• ■ • • - 

OUTP 

1 8 



OUTP 

• 19 — 



OUTP 

?.n 



OUTP 

21 * 


... « — 

OUTP 

22 



OUTP 

2 3 



OUTP 

24 



OUTP 

2 5 



OUTP 

26 



OUTP 

27 ' 



OUTP 

28 



OUTP 

29 


- 

OUTP 

30 




ICO FORMAT ( 3X* 3< I 2. IX) * 8X*F9*3. 1 0 X * A6 » 1 2 X * 3 A8 ) 
1C? FORMAT ( 20 X • F9# 3 • 1 0 X • A6 * t 2X * 5A 8 ) 

END 


OUTP— 31 

OUTP 32 
OUTP 33 
OUTP 34 



I. 


IDENTIFICATION 


A. Name: PRTOUT 

B. Language: FORTRAN IV 

C. Machine: IBM 360 


D. Purpose: To arrange the predicted antenna-noise 

temperature results for print and plot. 



E. Calling 

Sequence : 

CALL PRTOUT (TARY, KARY, NSTA, 
PLT, Kl) • 

Name 

Dimension 

T yP e - 

Description 

TARY 

(336,7) 

R*4 

Input-Array stored- the predicted 
temperatures of each stations. 

KARY 

(336,7) 

1*2 

Input-Indicate the temperature 
sources . 

NSTA 

CD 

1*4 

Input-Number of tracking ground 
stations . 

PLT 

CD 

L*1 

Input-"TRUE" , plot the results. 
"FALSE", no plot output. 

K1 

(D 

1*4 . 

Input-Tracking frequency within 
a tracking period. 


F. Common 

Block 


Name 

Variable 

Dimension 

Type Description 

/STAID/ 

NAME 

(7) 

R*8 Input-Name of the tracking 

stations . 


TYPE 

(7) 

1*4 Input-Tracking system type 

code . 

/MISCEL/ EPOCHY 

(D 

1*4 Input -Date (year, month, day) 

of orbit's epoch 
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Variable 

Dimension 

Type 

Des cription 

- - EPOCH! I 

(D 

1*4 

Input-Date (hour, minute, 
second) of orbit's epoch. 

MONA . 

(12) 

1*4 

Input-Name of each month, 




(4 characters) 

.MONT 

‘ (12) 

1*4 

Input-Cummulative days of 
a year. 


G. Non- System Routines Required: RYMDI , OUTPUT, TPLOT 
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LEVEL 13 < SEPT 6<J » 


OS/360 FORTRAN H 


DATE T2. 061/14, 44. SO 


MAIN# OP T=0 1 .L X NEC NT =58 , SOURCE .EBCDIC *NOL 1ST. NO DECK 


— - 

ISN 

0002 


ION 

o : c 3 


ISN 

0 C 0 4 


I SN 

Sue- *5 

- 

I S N 

noC6 


ISN 

ISM 

ni, j 7 

r: s •; a -- • 


— 

t SN 

e-;c9 


I jN 

c c : > 


ISM 

"•Ml 

• 

tr.N 

0012. 


ISN 

0 0 11 



ISM 

CC14 

- 

ISN 

0 0 1 s 


I S N 

5 0 1 6 


ISN 

s •: l r 

- 

r sn 

r " H 


ISN 

OC 19 



r sn 

0020 • — 

- 

ISN 

0 0 21 

— 

ISN 

C r 22 


I SN 

^ 23 

— — — 

ISN 

OC 24 


ISN 

ISN 

SC 25 


C 0 26 -• 

- 

ISN 

cr.p 7 

Cn 

ISN 

0 C 2 9 

O' 

ISN 

C C 30 

* 

ISN 

00 31 


ISN 

0 0 3 3 

— 

t SN 


- 

ISN 

or 35 



ISN 

CO 36 


t SN 

0 C 3 7 

3 

ISN 

0 0 33 


ISN 

OS 39 

— 

ISN 

CC4C • 

0 

l SN 

C 0 42 


ISN 

09 4 3 


1 SN 

(T4A 

vj 

t SN 

C04S 

ISN 

C C 4 7 


t SN 



I SN 

CO *9 ■ 

- ••• “■ 

ism 

0 0 50 


(SN 

C " 5 1 . 


t SN 

CO 52 


tSN 

tSN 

C C 54 


00 55 


ISN 

0 C 56 


tSM 

C 6 57 


2SN 

0 0 58 


COMPILER OPTIONS - NAME= 

C 

• - SUBROUTINE PRTOUT I TA9Y.KARY.NSTA.PLT.K1 ) ~ 

INTEGER*2 KC.KARYI 336.7) 

INTEGERT4 EPOCHY, EPOCHH. INTRVL. TYPE! T1 .DAY 
REAL*4 TAPYI 336.7) 

RCAL*8 NAMEI7) 

LOGICAL* 1 T I ME t PL T 

COMMON /STA 10/NAME , TYPE 

COMMON/MISCEL/EPOCHY .EPOCHH. INTRVL. ISTRTY.ISTRTH 
common /ALMN/M0N4I I 2 ) .MONTI 1 2) - ■ 

T I M E= .FALSE. 

CALL R YMO I I I STRTY , IR Y . IRM, I RO) - 
CALL RYMOKISTRTH.IHE.IME.ISEI 

• ! L- IRM • ■ 

DAY = IPO 
LOLO=IPM 
IOUAY= IPO 

NU’-'O Y = OA Y +MOnV I L“ l ) 

L INE= l 

WRITE (6. ICAO) MONAIL I .0AY.1RY 

WRITE I 6. 10 SO) 

I SEC= I HE *3 600+ I ME » 60 *1 SE- INTRVL 
•DO 30 I=1.K1 

1 5 E C = I SEC ♦ INTRVL 
CO 20 J-l.NSTA 

TA= TAPYI t , J) 

I E I TA.FO.O) GO TO 20 
K C = K A R Y ( t, J ) 

CALL OUt PUT I I SEC, NAME! J) .TA.KC.TIME) 

IFtLINE.GT.50) GO TO 21 
L I NL-L I NE* i 

00 to 2- 

21 WR!TE(6,1C4C> MONAIL). OAY.IRY 

VP 1 TEt6. 1 ;s?) 

■t IME-.TPUE* 

LINE=1 

20 CONTINUE 

IEI ISFCMNTRVL .LT. 66400) GO TO 3C ” 

I SEC= I SEC-064CO 
NUMDY=NUMOVF l 

fF( NUMDY.LE .MONT(L) > GO to 31 
L -L * 1 

31 0AY=NUMOY-MONTtL-l ) 

WRt TEt 6. 1040 ) MONAI L ) .OA Y, t RY ' 

-■* WR I TCI 6. t 0 53 ) 

L I N EM t 
30 COntInUE 

IF I .Nrtt.pLT) Return 

CALL tPLOTlLOLD.LitOOAV.bAY* t NT RVL. NSTAj TARY ) 

loot RdRMA+t tHl , 40X, A4.3X. t 2. •« * 4 3X« MR* *12) - -•••••- 

iCSC FORMAT t //, 5K« • T tME • , l0x« ‘TEMPERATURE • ,8* ♦ • TRACK 1 NG STATION* . 
i 5X4 ‘COMMENT* , /3 X 4 * HH MM SS ‘ . 7X i * KELV I N DEGREES • ) 

RETURN 

ENo > - 


10.NOXREF 



<3 


e> 





I. 


IDENTIFICATION 


A. 

Name : 

TPLOT 


B. 

Language 

: FORTRAN 

IV 

C. 

Machine : 

IBM 360 


D. 

Purpose : 

To plot the predicted Antenna-noise 
temperature on the Calcomp plotter. 

E. 

Calling 

Sequence : 

CALL PLOT (LOLD, LNEW , IODAY, 
I DAY, INTRVL, N, TARY) 


Name 

Dimension 

Type 

Description 

LOLD 

(1) 

I 

Input-Month of the starting date 
of predicted period. 

LNEW 

(« 

I 

Input-Month of the end date of 
predicted period. 

IODAY 

(D 

I 

Input-Day of the starting date. 

I DAY 

(D 

I 

Input-Day of the end date. 

INTRVL 

(1) 

I 

Input-Tracking interval (seconds) 

N 

(1) 

I 

Input-Number of the tracking 
stations. 

TARY 

CD 

I 

Input-Predicted temperature data 
array to be plot. 


F . Common 

Blocks : 


Name 

Variable Dimension 

Type ' Description 

/Aim/ 

MONA 

12 

I Input-Name of each month 


MONT 

12 

I Input-Cummulative days of 

each month 

/STAID/ 

NAME 

m 

R*8 Input-Tracking station 

identification 
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Name 


Variable Dimension Type Description 



TYPE 

(7) 

1*4 

Input-Tracking antenna 
type code . 


LATD 

(7) 

1*4 ) 

Input -Degrees , Minutes . 


LATM 

(7) 

1*4 > 

and seconds of station's 


LSLAT 

(7) 

R*4 ) 

geodetic latitude. 


LOND 

(7) 

1*4 ) 

Input Degreees , Minutes 


LONM 

(7) 

1*4 > 

and Seconds of station's 


LSLONG 

(7) 

1*4 ) 

geodetic longitude. 


LGAINA 

(7) 

R*4 

Input-Receiver’s peak 
antenna gain (dB) 


LELMIN 

(1) 

R*4 

Input-Minimum observation 
altitude angle (deg.) 


NSTA 

(1) 

1*4 

Input-Number of tracking 
stations 

/ANTNA/ 

ANTEN 

(3) 

to 

X 

oo 

Tracking antenna type. 

/NEW/ 

FREQ 

(1) 

1*4 

Transmitter frequency 
(megacycles /second) 


Non-System Routines 

Required : 

Name 

Entry Point • 

CALCOMP 

FRMADV 

EDIT 

EDIT 

GRID 

GRID, SETGRD , PLOT 

HORLIN 

HORLIN, VERLIN 

PLOST 

PLOST, ENDPLT 


III. METHOD 

Subroutine TPLOT utilizes a plot package developed by 
WOLF Research and Development Corporation. Hardware options 
for plotted output include the CALCOMP plotter or computer 
printer, or both. 
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The CALCOMP plotter is specified in the main program 
and plotted output with associated plotter commands are stored 
on a 7- track magnetic tape, mounted on logical unit PLOTAPE. 

Predicted telemetry data were stored in TARY . The 
predicted time array corresponding to predicted telemetry 
data were stored in XARY . 

The following operations are performed by subroutine 
TARY in producing predicted data plots: 

1. Conversion of prediction time (seconds) of each 
plotted value to the zero hour reference frame 
of the first plotted point, when a pass extends 
beyond a day. 

2. Determination of the density of vertical grid 
lines from the time length of predicted period. 

3. Conversion of time to days. 

4. Determination of the density of horizontal grid 
lines from the maximum and minimum plotted value. 

5. Identification of plot frame. 

6. Plotting of antenna-noise temperature. 
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LEVEL IS '< SEPT 69 £1 


'OS/369 FORTRAN H 



TSN 

SCO 2 


ISN 

9 CC 3 


TGN 

0 0 0 4 


T5N 

C 0 C 5 

n 

•ISN 

">2 06 


ISN 

0 00 7 



ISN 

CO j* 

#% 

' • ' ■ - • ISN 

? r S 9 


r sn 

00 1 0 


' J SN 

OCl '2 

<s> 

ISN 

OC 12 


ISN 

C 0 l 3 


- ' :iSN 

ri Cl 4 


ISN 

Z CIS 


ISn 

00 16 

** 

ISN 

0 017 


ISN 

00 18 


;isn 

0 C 20 

*> 

• ;• isn 

CO 21 


T SN 

0)22 


.j ,, N 

C 0 2 3 


ISN 

CO 24 

•*-0 

... _ .. .j SN . 

:o as 

o 

I SN 

0026 


ISN 

0 C 27 


'ISN 

00 2 ^ 


— ■ ‘ 1 £N 

C r 29 

-- 

'ISN 

0 J T ) 


‘ISN 

on. ii 


ISN 

0 **• 32 


' • ISN 

00 33 


ISN 

''C 3-1 


ISN 

o-: 3 s 


;i r.N 

0 0 36 


I SN 

0 0 37 


I SN 

r. 0 3 B 


ISN 

0 .: 39 


I S N 

0 0 40 


ISN 

0 J 4 'l 

S 

I SN 

? J '42 


‘ ISN 

0043 


ISN 

CO 44 

S' 

ISN 

0 045 


ISN 

OC 46 


ISN 

C 0 4 7 

f • 

ISN 

0V»8 


*"* ' ISN 

0*' 4 9 


, ISN 

0 0 50 


ISN 

or 51 


^PILER OPTTONS — NYAME-= SHA'XN*# OPT=0*lV LI NEC NT = 5 8 • SOURCE * E8CD I C • 
C 

... SUBROUTINE TPLOTC LOLO »LN£ W * I OD AY • IOAY*lNTRVL«N*T ARY ) 

DIMENSION MQNf( 12). XAPYU36) 

:tNTCGEW*4 TNTG'4 * T YPF ( 7 ) ,FngQ 

;PE- ALVA LSL'A r< T) ;lSLON( 7) ,TARY< 336.7 ) #LELMI N.LGA X NA 
REALM'S NAME < 7) -,'OUMY, W0RD3 . A NTE N< 3 > 

LOGTCAL*’l mr-TLi:.( 36) /• 1973 TO , 

* TirTL2( 339'/ 'TRACKING ’INTERVAL": MINUTES*'/ * 

■» T;!TL3( S'!') /’•‘TRACKING STATION ;t ANTENNA 

A TTTL 4<1 9) / 'FREOUENCY 5 M HZ*/ 

LOG ICAL*’! BYTES(8) .OYT( A > ,0UT(4 ) *PCH 
' O AT A .2*3. 31.39*31* 39 .31 *31*30*31 *3 0*31/ 

EOUT VALENCE ( W0RD8, BYTES ( 1’))-. ( I NTG4 , QYT{ !)) 

COMMON 7ALMN/MONA ( 12) .MONTHS) 

COMMON'/ ST'A'I D/NAME »TYP£, L ATD < 71 ;LATM(7) * L SL AT . LONO< 7 » 
S. L SLON#LGATNA( 7) *t.ELMI N.NSTA 
COMMON //ANTNA/ANTEN 
COMMON /NEW/F R£Q«PCH 
T)AYR-IOOAY 

OAYP = OAYR*7 • - : 

TMA = aCCC « 0 

IF(:FRE0«E0*-400'I TMA = 20C0»C 
'I NT» = TMA/4, 

■X‘1 kFLO’ATOI NTRVL'I /064OO# 

JK = 7-./Xl 

: DO n-4 <=:i-,UK : ;*** 

14 X AR Y.(K )sf(X“T )*X lM-OAYR 

DO 12 J=’1-*N . ..... 

•DO >15 l=:i • JK 

‘VARY (T-» 'J ) - A M I N.l'( T ARY (T-» 'J j * TMA) 

15 CONTINUE 

'CALL FPMAOV - ’ 

CALL SFiT0f>D*t:UvlWl0**8*,-4) 

•CALL GRID (0AYRV0AYP.7** ‘•-*0 • 0**tTMA *4V* •-.0 tO) ‘ * 

CALL SCALE! t)AYR ,t>AYP. 0* * TMA ,0) 

CALL HOPLINI 34HRAE-B TRACKING ANTE NN A TE MPER ATURE • 34 
I N f G4= MOV A ( LOLD ) 

2C t i tl r( t. ) -r y t (u*) 

I Q A YR-'f F'i-X ( ! .OA Yfl*) 

C ALL ED I T ( I OAYR • •121 • *T I TLT-C 7) .NN.IBL'l 

INT CA=MO.NA'<'LNEW 1 
00 22 ’ L 2 ~-l**"4 

22 T I'TLTt 2 1*L 2 ) =GYT( L2 ) ' 

NUMBER- IOAY-T 

C ALL EO I t ( NUMBER.’* 12)*# 71 TL'1( 2 8 ) • NN* I BL ) 

NUMSER='INTP VL/60 

CALL EO I T ( NUMBER • *14) • t T I TL2 < 22 > *NN • f BL > 

WORDQ = NAME< J ) 

DO 40 I J = 1 , 6 ~r — - - 

4C T I TL3( 20M J )=BYTES< I J) 

WGRD0 = AN1cN( TYPE! J) ) ^ — "~7 

on 4i ij=4,b I 

41 TITL3I43*IJ|=BYTES(TJ) to 


<OATE 7.5? •061/14*45*00 ' 

NOL 1ST* NODECK* LOAD . MAP • NOED IT * ID. NO XREF -- • - 
TPLOT 

TPLO 1— 


• LONM ( 7 ) » 





TPLO 
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TPLO 
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6 


TPLO- 

- 7 * 

e 

TPLO 

e 


TPLO 

Q 


TPLO 

10 


TPLO 

1 1 


TPLO 

12 

— 

TPLO- 

13- 


TPLO 

14 


TPLO 

15 


TPLO 

1 6 


TPLO 

1 7 


TPLO 

1 8 


TPLO ' 

19 - 


TPLO 

2"* 


TPLO 

21 


TPLO 

22 


TPLO 

23 


TPLO 

24 

— 

TPLO- 

25 


TPLO 

26 


TPLO 

27 


TPLO 

28 


TPLO 

29 


TPLO 

30 


TPLO 

31 


TPLO 

32 


TPLO 

3.1 


TPLO 

34 

) 

t pi n 

35 


TPLO 

36 


TPLO 

37 


TPl.O 

38 


TPLO 

39 


TPLO 

40 


TPLO 

41 


TPLO 

42 
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ISN 3052 
ISN 0 C 53 

— ISN C95 4 

ISN GOSS 
ISN 0>55 
ISN * 0 57 

ISN 0 ^ 58 
ISN 3C59 

ISN COS'* 

ISN ->r.6l 

- — — ■ ISN ? 062 

ISN *C53 
ISN ;064 
ISN 0066 

isn CC63 

ISN CJ69 
— ‘ • ISN S3 73 

ISN O', 7 1 
ISN 0072 
ISN 0C73 

ISN 3074 

ISN 0 C 7 6 
ISN 0077 
ISN 3 -i 73 

ISN 0070 

ISN 0360 
; ISN 30 31 


CAUL E0tT( FREQ* *1 3) * .TI TL4< 1 3) ,NN,t BL) 

CALL HORL1 N( T ITL1 * 36,5. 0 *8* 3 .0 *0 ) 

CALL HORLI N( TI TL 2 . 33 . 3« 5 .8. I .0 *0 ) -*•” 

CALL HORL IN( T I TL 3 « 5 1 « 5* O • 8* 5 • 0 • O ) 

CALL HOPLIN(TI TL4, 1 9 • 6* 5 . 0. 1 . 0 . C > 

CALL VFPLINt'T E M P E R A TUR E * • 2 1 « 0 • 1 • 4 o 5 * 0 *0 ) 
CALL VERLINI •KELVIN DEG • • * 1 1 ,0 • 3 * 4 • 5 . C « 0 > 

CALL HOPLtN(«TtME (DAY OF THE MONTH ) • « 23 • 5 • 0 « 0 « 7 *0 • 0 > 

K»n*C- C N1 - 1 > + 1*75 

NUM3EP=ZFI X< TMA >-( Nl- 1 ) *INTR 

CALL EO IT ( NUMBER . M 5> * .OUT.NN*IOLl 

|F( NUM0ER.E0.530) IRL^IBL-l 

IF( NUMOER*EO.O ) IUL=I0L-3 

60 CALL HO PL l N C OUT ( I 8L ) *NN*0#65*X*0*01 

DO 61 N 2=1*6 
X=1 «-(N2-l ) + 1*2857 
NUMBER "O AYR ♦ < N2- 1 ) 

NUMBER sMDOC NUMBER- 1 , MQN (LOLO) >4-1 
call EO IT (NUMBER* M 3) • • OUT • NN, I BL ) 

IF( NUMOER.LT* 10) I8L=I8L-t ** ~ 

6 t CALL HORLIN(OUT( I OL ) *NN * X,0* 9*0 .0 I 

CALL PLOT < XARY.TARY< 1 » J) • JK * • • ) ' - ~ 

12 CONTINUE 
CALL PAWS 
RETURN 
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